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Nomenclature 

A fuselage cross-sectional area 

Ag fuselage surface area 

A p frame cross-sectional area 

(. AR ) aspect ratio of wing 

b wingspan; intercept of regression line 

b s stiffener spacing 

bs wing structural semispan, measured along 

quarter chord from fuselage 

b w stiffener depth 

Cp Shanley’s constant 

Cp center of pressure 

Cp root chord of wing at fuselage intersection 

Cp theoretical root chord of wing 

C s j portion of wing leading edge not used for 

structural box 

C Sl portion of wing trailing edge not used for 

structural box 

Csp structural root chord of wing 

Cst structural tip chord of wing 

Cp tip chord of wing 

d frame spacing 

d\\/ optimum web spacing of wing 

D maximum diameter of fuselage 

e wing buckling exponent 

ec wing cover material factor 

E Young’s modulus of shell material 

Ep Young’s modulus of frame material 

F cv compressive yield strength 

Fs shear strength 

F tu ultimate tensile strength 

h thickness of sandwich shell 

h e . step function for i th engine on wing 

hfgj step function for i th landing gear on wing 

Ip frame cross-sectional area moment of inertia 

7 V area moment of inertia about the y-axis 

/y V 


2 

Kpj frame stiffness coefficient, IpIAp 

K m g shell minimum gage factor 

Kp shell geometry factor for hoop stress 

K$ constant for shear stress in wing 

K t h sandwich thickness parameter 

ip fuselage length 

iLE length from leading edge to structural box at 

theoretical root chord 

l MG length from nose to fuselage mounted main 

gear 

/yVG length from nose to nose gear 

Ipp length from trailing edge to structural box at 

theoretical root chord 

/ ] length of nose portion of fuselage 

12 length of tail portion of fuselage 

l n length from nose to breakpoint of fuselage 

L lift 

Lp maximum vertical tail lift 

m buckling equation exponent; slope of 

regression line 

M longitudinal bending moment 

n normal load factor 

n x longitudinal acceleration 

N Xa axial stress resultant 

N Xb bending stress resultant 

N xp pressure stress resultant 

N x tensile axial stress resultant 

N~ compressive axial stress resultant 

N y hoop direction stress resultant 

P perimeter 

Pg internal gage pressure 

P s perimeter of shell 

P w perimeter of walls 

P | exponent of power law of nose section of 

fuselage 

P 2 exponent of power law of tail section of 

fuselage 

r radius of fuselage 



r(y) total wing chord as a function of position 

along quarter chord 

r s (y 0 structural wing chord as a function of position 

along quarter chord 

R correlation coefficient used for regression 

Rjj n fineness ratio 

Rht ratio of horizontal tail station to fuselage 
length 

RlE ratio °f wing leading edge station at 

theoretical root chord to fuselage length 

RmG ratio °f length to main gear to fuselage length 

Rnq ratio of length to nose gear to fuselage length 

Rp j ratio of length to leading edge of fuselage 

mounted propulsion to fuselage length 

Rp2 ratio of length to trailing edge of fuselage 

mounted propulsion to fuselage length 

Rtiy) thickness ratio of wing as a function of 
position along quarter chord 

RjAP taper ratio of wing 

Sp plan area of the fuselage 

Spc stroke of landing gear 

Sp plan area of wing 

t( y) thickness of wing box as a function of position 

along quarter chord 

t c core thickness 

tj face sheet thickness 

tg material gage thickness, t s I K m g 

t m g material minimum gage thickness 

t s skin thickness 

t w stiffener thickness 

t total equivalent isotropic thickness of shell and 

frames 

Ip total equivalent isotropic thickness of fuselage 

structure 

tp smeared equivalent isotropic thickness of 

frames 


h c 

shell thickness required to preclude 
compressive failure 

h G 

shell thickness required to meet minimum 
gage constraint 

hr 

shell thickness required to preclude tensile 
failure 

h 

smeared tension tie thickness 

i\v 

smeared wall thickness 

he 

thickness of wall to meet minimum gage 
constraint 

t\Vj 

thickness of wall required to prevent tensile 
failure 

T 

torque on wing carrythrough structure 

VB 

fuselage volume 

Vw 

volume of wing structural box, including 
structural components 

V\ 

volume of nose section of fuselage 

v 2 

volume of tail section of fuselage 

W C 

width of carrythrough structure of wing 

w 

weight of aircraft structure 

W' 

weight of wing per unit span 

W B 

weight of fuselage structure and attached 
components 

Wpp 

weight of fuel 

Wi 

ideal fuselage structural weight 

Wno 

weight of nonoptimum material 


vehicle longitudinal weight distribution 

Wto 

gross takeoff weight of aircraft 

w/s 

shell structural weight per unit surface area 

X 

longitudinal fuselage coordinate 

x calc 

weight calculated by PDCYL 

X HT 

distance from nose to theoretical quarter chord 
of horizontal tail 

X LE 

distance from nose to leading edge of wing at 
theoretical root chord 


f S 

Isp 


equivalent isotropic thickness of shell 

shell thickness required to preclude buckling 
failure 


X P i 


xp 2 


distance from nose to leading edge of fuselage 
mounted propulsion 

distance from nose to trailing edge of fuselage 
mounted propulsion 


IV 



y 

transverse fuselage coordinate; wing 

zw 


coordinate measured along quarter chord 

A 

}'act 

actual weight 

p 

yest 

estimated weight after regression 

p 

V 

vertical fuselage coordinate 

r 

Z(y) 

total width of wing box as a function of 

r I) 


position along quarter chord 

PF 

Zs(y) 

width of wing box structure as a function of 

OS 


position along quarter chord 

1 

8 

frame deflection 


E 

shell buckling efficiency 

V 

ec 

wing cover structural efficiency 



wing web structural efficiency 
wing sweep 
wing loading 

structural material density 

gross fuselage density 

frame structural material density 

allowable shear stress for wing 

sum over fuselage or wing length; solidity of 
wing 

truss core angle 


v 
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Summary 

A method of estimating the load-bearing fuselage weight 
and wing weight of transport aircraft based on funda- 
mental structural principles has been developed. This 
method of weight estimation represents a compromise 
between the rapid assessment of component weight using 
empirical methods based on actual weights of existing 
aircraft, and detailed, but time-consuming, analysis using 
the finite element method. The method was applied to 
eight existing subsonic transports for validation and corre- 
lation. Integration of the resulting computer program, 
PDCYL, has been made into the weights-calculating 
module of the AirCraft SYNThesis (ACSYNT) computer 
program. ACSYNT has traditionally used only empirical 
weight estimation methods; PDCYL adds to ACSYNT a 
rapid, accurate means of assessing the fuselage and wing 
weights of unconventional aircraft. PDCYL also allows 
flexibility in the choice of structural concept, as well as a 
direct means of determining the impact of advanced 
materials on structural weight. 

Using statistical analysis techniques, relations between 
the load-bearing fuselage and wing weights calculated by 
PDCYL and corresponding actual weights were deter- 
mined. A User’s Manual and two sample outputs, one for 
a typical transport and another for an advanced concept 
vehicle, are given in the appendices. 

Introduction 

A methodology based on fundamental structural 
principles has been developed to estimate the load- 
carrying weight of the fuselage and basic box weight of 
the wing for aircraft, and has been incorporated into the 
AirCraft SYNThesis program (ACSYNT). This weight 
routine is also available to run independently of 
ACSYNT, and is a modification of a collection of pre- 
viously developed structural programs (refs. 1-4). The 
main subroutine called by ACSYNT is PDCYL. This 
study has concentrated on modern transport aircraft 
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because of the detailed weight information available, 
allowing the weights output from PDCYL to be compared 
to actual structural weights. The detailed weight state- 
ments also allow nonoptimum factors to be computed 
which, when multiplied by the load-bearing structural 
weights calculated by PDCYL, will give good representa- 
tive total structure weight estimates. These nonoptimum 
factors will be computed through a regression analysis of 
a group of eight transport aircraft. 

PDCYL is able to model both skin-stringer-frame and 
composite sandwich shell fuselage and wing box 
constructions. Numerous modifications were made to 
PDCYL and its associated collection of subroutines. 

These modifications include the addition of detailed 
fuselage shell geometry calculations; optional integration 
of a cylindrical fuselage midsection between the nose and 
tail sections; addition of landing and bump maneuvers to 
the load cases sizing the fuselage; ability to introduce an 
elliptical spanwise lift load distribution on the wing; 
variation of wing thickness ratio from tip to root; ability 
to place landing gear on the wing to relieve spanwise 
bending loads; distribution of propulsion system compo- 
nents between wing and fuselage; and the determination 
of maximum wingtip deflection. 

Brief Description of ACSYNT 

The Aircraft Synthesis Computer program, ACSYNT, 
is an integrated design tool used in the modeling of 
advanced aircraft for conceptual design studies (ref. 5). 
ACSYNT development began at NASA Ames Research 
Center in the 1970s and continues to this day. The 
ACSYNT program is quite flexible and can model a wide 
range of aircraft configurations and sizes, from remotely 
piloted high altitude craft to the largest transport. 

The ACSYNT program uses the following modules, not 
necessarily in this order: Geometry, Trajectory, Aero- 
dynamics, Propulsion, Stability, Weights, Cost, Advanced 
Aerodynamic Methods, and Takeoff. An ACSYNT run 
would normally progress as follows: the Geometry 
module is called to define the aircraft shape and configu- 
ration; the Trajectory module then runs the vehicle 
through a specified mission; finally the Weight and Cost 



modules are executed. To determine the performance of 
the vehicle at each mission point, the Trajectory module 
will call the Aerodynamics and Propulsion modules. 

After the mission is completed, the calculated weight of 
the aircraft may be compared with the initial estimate and 
an iteration scheme run to converge upon the required 
aircraft weight. This process is necessarily iterative as the 
aircraft weight ACSYNT calculates is dependent upon the 
initial weight estimate. 

ACSYNT is able to perform a sensitivity analysis on any 
design variable, such as aspect ratio, thickness-to-chord 
ratio, fuselage length or maximum fuselage diameter. 
Sensitivity is defined as (change in objective function/ 
value of objective function) divided by (change in design 
variable/design variable). As an example, if gross weight 
is the objective function and decreases when the wing 
thickness-to-chord ratio increases, then the sensitivity of 
thickness-to-chord ratio is negative. It is important to note 
that while this increase in thickness-to-chord ratio lowers 
the gross weight of the aircraft, it may also have a 
detrimental effect on aircraft performance. 

ACSYNT is also able to size multiple design variables by 
optimizing the objective function. The objective function 
represents the interactions between design disciplines 
such as structures, aerodynamics and propulsion. The 
automated sizing of design variables during the optimi- 
zation process is accomplished using the gradient method. 
Two types of constraints may be imposed during the 
optimization process. These are performance-based 
constraints such as runway length or maximum roll angle, 
and side constraints on design variables such as limita- 
tions on wing span or fuselage length. ACSYNT never 
violates constraints during the optimization process so 
that each iteration produces a valid aircraft. 

Methods of Weight Estimation 

Two methods are commonly available to estimate the 
load-bearing fuselage weight and wing box structure 
weight of aircraft. These methods, in increasing order of 
complexity and accuracy, are empirical regression and 
detailed finite element structural analysis. Each method 
has particular advantages and limitations which will be 
briefly discussed in the following sections. There is an 
additional method based on classical plate theory (CPT) 
which may be used to estimate the weight of the wing box 
structure. 

Empirical- The empirical approach is the simplest 
weight estimation tool. It requires knowledge of fuselage 
and wing weights from a number of similar existing 
aircraft in addition to various key configuration parame- 
ters of these aircraft in order to produce a linear regres- 


sion. This regression is a function of the configuration 
parameters of the existing aircraft and is then scaled to 
give an estimate of fuselage and wing weights for an 
aircraft under investigation. Obviously, the accuracy of 
this method is dependent upon the quality and quantity of 
data available for existing aircraft. Also, the accuracy of 
the estimation will depend on how closely the existing 
aircraft match the configuration and weight of the aircraft 
under investigation. All of the empirical regression 
functions currently in the ACSYNT program give total 
fuselage weight and total wing weight. 

Finite Element- Finite element analysis is the matrix 
method of solution of a discretized model of a structure. 
This structure, such as an aircraft fuselage or wing, is 
modeled as a system of elements connected to adjacent 
elements at nodal points. An element is a discrete (or 
finite) structure that has a certain geometric makeup and 
set of physical characteristics. A nodal force acts at each 
nodal point, which is capable of displacement. A set of 
mathematical equations may be written for each element 
relating its nodal displacements to the corresponding 
nodal forces. For skeletal structures, such as those 
composed of rods or beams, the determination of element 
sizing and corresponding nodal positioning is relatively 
straightforward. Placement of nodal points on these 
simple structures would naturally fall on positions of 
concentrated external force application or joints, where 
discontinuities in local displacement occur. 

Continuum structures, such as an aircraft fuselage or 
wing, which would use some combination of solid, flat 
plate, or shell elements, are not as easily discretizable. An 
approximate mesh of elements must be made to model 
these structures. In effect, an idealized model of the 
structure is made, where the element selection and sizing 
is tailored to local loading and stress conditions. 

The assembly of elements representing the entire structure 
is a large set of simultaneous equations that, when com- 
bined with the loading condition and physical constraints, 
can be solved to find the unknown nodal forces and 
displacements. The nodal forces and displacements are 
then substituted back into the each element to produce 
stress and strain distributions for the entire structural 
model. 

Classical Plate Theory- CPT has been applied to wing 
structure design and weight estimation for the past 
20 years. Using CPT a mathematical model of the wing 
based on an equivalent plate representation is combined 
with global Ritz analysis techniques to study the struc- 
tural response of the wing. An equivalent plate model 
does not require detailed structural design data as required 
for finite element analysis model generation and has been 
shown to be a reliable model for low aspect ratio fighter 
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wings. Generally, CPT will overestimate the stiffness of 
more flexible, higher aspect ratio wings, such as those 
employed on modern transport aircraft. Recently, 
transverse shear deformation has been included in 
equivalent plate models to account for this added 
flexibility. This new technique has been shown to give 
closer representations of tip deflection and natural 
frequencies of higher aspect ratio wings, although it still 
overestimates the wing stiffness. No fuselage weight 
estimation technique which corresponds to the equivalent 
plate model for wing structures is available. 

Need for Better, Intermediate Method 

Preliminary weight estimates of aircraft are traditionally 
made using empirical methods based on the weights of 
existing aircraft, as has been described. These methods, 
however, are undesirable for studies of unconventional 
aircraft concepts for two reasons. First, since the weight 
estimating formulas are based on existing aircraft, their 
application to unconventional configurations (i.e., canard 
aircraft or area ruled bodies) is suspect. Second, they 
provide no straightforward method to assess the impact 
of advanced technologies and materials (i.e., bonded 
construction and advanced composite laminates). 

On the other hand, finite-element based methods of 
structural analysis, commonly used in aircraft detailed 
design, are not appropriate for conceptual design, as the 
idealized structural model must be built off-line. The 
solution of even a moderately complex model is also 
computationally intensive and will become a bottleneck 
in the vehicle synthesis. Two approaches which may 
simplify finite-element structural analysis also have draw- 
backs. The first approach is to create detailed analyses at 
a few critical locations on the fuselage and wing, then 
extrapolate the results to the entire aircraft, but this can be 
misleading because of the great variety of structural, load, 
and geometric characteristics in a typical design. The 
second method is to create an extremely coarse model of 
the aircraft, but this scheme may miss key loading and 
stress concentrations in addition to suffering from the 
problems associated with a number of detailed analyses. 

The fuselage and wing structural weight estimation 
method employed in PDCYL is based on another 
approach, beam theory structural analysis. This results 
in a weight estimate that is directly driven by material 
properties, load conditions, and vehicle size and shape, 
and is not confined to an existing data base. Since the 
analysis is done station-by-station along the vehicle 
longitudinal axis, and along the wing structural chord, the 
distribution of loads and vehicle geometry is accounted 
for, giving an integrated weight that accounts for local 
conditions. An analysis based solely on fundamental 


principles will give an accurate estimate of structural 
weight only. Weights for fuselage and wing secondary 
structure, including control surfaces and leading and 
trailing edges, and some items from the primary structure, 
such as doublers, cutouts, and fasteners, must be esti- 
mated from correlation to existing aircraft. 

The equivalent plate representation, which is unable to 
model the fuselage structure, is not used in PDCYL. 

Methods 

Overview 

Since it is necessary in systems analysis studies to be able 
to rapidly evaluate a large number of specific designs, the 
methods employed in PDCYL are based on idealized 
vehicle models and simplified structural analysis. The 
analyses of the fuselage and wing structures are per- 
formed in different routines within PDCYL, and, as such, 
will be discussed separately. The PDCYL weight analysis 
program is initiated at the point where ACSYNT per- 
forms its fuselage weight calculation. PDCYL first 
performs a basic geometrical sizing of the aircraft in 
which the overall dimensions of the aircraft are deter- 
mined and the propulsion system, landing gear, wing, and 
lifting surfaces are placed. 

Fuselage- The detailed fuselage analysis starts with a 
calculation of vehicle loads on a station-by-station basis. 
Three types of loads are considered — longitudinal 
acceleration (applicable to high-thrust propulsion 
systems), tank or internal cabin pressure, and longitudinal 
bending moment. All of these loads occur simultaneously, 
representing a critical loading condition. For longitudinal 
acceleration, longitudinal stress resultants caused by 
acceleration are computed as a function of longitudinal 
fuselage station; these stress resultants are compressive 
ahead of the propulsion system and tensile behind the 
propulsion system. For internal pressure loads, the 
longitudinal distribution of longitudinal and circumferen- 
tial (hoop) stress resultants is computed for a given shell 
gage pressure (generally 12 psig). There is an option to 
either use the pressure loads to reduce the compressive 
loads from other sources or not to do this; in either case, 
the pressure loads are added to the other tensile loads. 

Longitudinal bending moment distributions from three 
load cases are examined for the fuselage. Loads on 
the fuselage are computed for a quasi-static pull-up 
maneuver, a landing maneuver, and travel over runway 
bumps. These three load cases occur at user-specified 
fractions of gross takeoff weight. Aerodynamic loads are 
computed as a constant fraction of fuselage plan form area 
and are considered negligible for subsonic transports. For 
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pitch control there is an option to use either elevators 
mounted on the horizontal tail (the conventional config- 
uration) or elevons mounted on the trailing edges of the 
wing. The envelope of maximum bending moments is 
computed for all three load cases and is then used to 
determine the net stress resultants at each fuselage station. 

After the net stress resultants are determined at each 
fuselage station, a search is conducted at each station to 
determine the amount of structural material required to 
preclude failure in the most critical condition at the most 
critical point on the shell circumference. This critical 
point is assumed to be the outermost fiber at each station. 
Failure modes considered are tensile yield, compressive 
yield, local buckling, and gross buckling of the entire 
structure. A minimum gage restriction is also imposed as 
a final criterion. It is assumed that the material near the 
neutral fiber of the fuselage (with respect to longitudinal 
bending loads) is sufficient to resist the shear and torsion 
loads transmitted through the fuselage. For the shear 
loads this is a good approximation as the fibers farthest 
from the neutral axis will carry no shear. Also, for beams 
with large fineness ratios (fuselage length/maximum 
diameter) bending becomes the predominant failure 
mode. 

The maximum stress failure theory is used for predicting 
yield failures. Buckling calculations assume stiffened 
shells behave as wide columns and sandwich shells 
behave as cylinders. The frames required for the stiffened 
shells are sized by the Shanley criterion. This criterion is 
based on the premise that, to a first-order approximation, 
the frames act as elastic supports for the wide column 
(ref. 6). 

There are a variety of structural geometries available for 
the fuselage. There is a simply stiffened shell concept 
using longitudinal frames. There are three concepts with 
Z-stiffened shells and longitudinal frames; one with 
structural material proportioned to give minimum weight 
in buckling, one with buckling efficiency compromised to 
give lighter weight in minimum gage, and one a buckling- 
pressure compromise. Similarly, there are three truss-core 
sandwich designs, two for minimal weight in buckling 
with and without frames, and one a buckling-minimum 
gage compromise. 

It is assumed that the structural materials exhibit elasto- 
plastic behavior. Further, to account for the effects of 
creep, fatigue, stress-corrosion, thermal cycling and 
thermal stresses, options are available to scale the 
material properties of strength and Young’s modulus 
of elasticity. In the numerical results of this study, all 
materials were considered elastic and the full room- 
temperature material properties were used. 


Composite materials can be modeled with PDCYL by 
assuming them to consist of orthotropic lamina formed 
into quasi-isotropic (two-dimensionally, or planar, 
isotropic) laminates. Each of the lamina is assumed to be 
composed of filaments placed unidirectionally in a matrix 
material. Such a laminate has been found to give very 
nearly minimum weight for typical aircraft structures. 

Wing- The wing structure is a multi-web box beam 
designed by spanwise bending and shear. The wing- 
fuselage carrythrough structure, defined by the wing- 
fuselage intersection, carries the spanwise bending, shear, 
and torsion loads introduced by the outboard portion of 
the wing. 

The load case used for the wing weight analysis is the 
quasi-static pull-up maneuver. The applied loads to the 
wing include the distributed lift and inertia forces, and the 
point loads of landing gear and propulsion, if placed on 
the wing. Fuel may also be stored in the wing, which will 
relieve bending loads during the pull-up maneuver. 

The wing weight analysis proceeds in a similar fashion to 
that of the fuselage. The weight of the structural box is 
determined by calculating the minimum amount of 
material required to satisfy static buckling and strength 
requirements at a series of spanwise stations. The covers 
of the multi-web box are sized by buckling due to local 
instability and the webs by flexure-induced crushing. 
Required shear material is computed independently of 
buckling material. Aeroelastic effects are not accounted 
for directly, although an approximation of the magnitude 
of the tip deflection during the pull-up maneuver is made. 
For the carrythrough structure, buckling, shear, and 
torsion material are computed independently and 
summed. 

As for the fuselage, there are a variety of structural 
geometries available. There are a total of six structural 
concepts, three with unstiffened covers and three with 
truss-stiffened covers. Both cover configurations use 
webs that are either Z-stiffened, unflanged, or trusses. 

Geometry 

Fuselage- The fuselage is assumed to be composed of a 
nose section, an optional cylindrical midsection, and a tail 
section. The gross density and fineness ratio are defined 
as 


P B = 

W B 

Vb 

(1) 

II 

l JL 

D 

(2) 
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where Wg is the fuselage weight (Wg = gross takeoff 
weight excluding the summed weight of the wing, tails, 
wing-mounted landing gear, wing-mounted propulsion, 
and fuel if stored in the wing), Vg is the total fuselage 
volume, Ig is the fuselage length, and D is the maximum 
fuselage diameter. The fuselage outline is defined by two 
power-law bodies of revolution placed back-to-back, with 
an optional cylindrical midsection between them (fig. 1). 
(For the present study, all eight transports used for 
validation of the analysis used the optional cylindrical 
midsection.) 

With the cylindrical midsection, integration gives the 
fuselage volume, fuselage planform area, and fuselage 
surface area as 


Vb 


kD 


2P\ + 1 


+(/ B -/ 2 -/|)+ 


l 2 


2Pl + \ 


(3) 


S B = D 


f\+ 1 


+ (/ B -/ 2 -/]) + 


h 


P2 + 1 


(4) 


A b = nS B 


(5) 


respectively, where l\ and / 2 are the respective lengths to 
the start and end of the cylindrical midsection, and P\ and 
P2 are the respective powers that describe the nose and 
tail sections. P\ and P 2 , again for the case of the cylin- 
drical midsection, arc found by solving the power-law 


equations for the volumes of the nose and tail sections, 
which are input from ACSYNT. The solution of these 
equations gives the respective nose and tail powers as 


_ 71 D 2 /, 1 

n 8V, 2 

(6) 

A - rcZ)2/ 2 > 

^ 8V 2 2 

(7) 


where V\ and V 2 are the corresponding nose and tail 
volumes. 

The horizontal tail is placed according to its quarter chord 
location as a fraction of the fuselage length. The distance 
from the nose to the tail is 

X HT = l B R HT (8) 

where Rj^j is the ratio of horizontal tail station to 
fuselage length. 

Propulsion may be either mounted on the fuselage or 
placed on the wing. In the case of fuselage mounted 
propulsion, the starting and ending positions of the 
propulsion unit are again calculated from their respective 
fractions of fuselage length as 


X P\ =l B R F\ 

(9) 

X P 1 ~ l B R P 2 

(10) 
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where Rp\ and Rp2 are the corresponding ratios of 
lengths to the leading and trailing edges of the fuselage 
engine pod to fuselage length. 

Similarly, the nose landing gear is placed on the fuselage 
as a fraction of vehicle length; the main gear, on the other 
hand, may be placed either on the fuselage as a single 
unit, also as a fraction of fuselage length, or on the wing 
in multiple units as will be described below. The positions 
of the respective nose and optional fuselage-mounted 
main gear are 

l NG- l B R NG 0 0 

l MG =l B R MG ( 12 ) 

where fijvc and R MG are the corresponding length ratios 
for the nose gear and main gear stations to vehicle length. 

Wing- The lifting planforms are assumed to be tapered, 
swept wings with straight leading and trailing edges. The 
planform shape is trapezoidal as the root chord and tip 
chord are parallel. 


The wing loading is defined as 


ft = 


W TO 

S P 


(13) 


where Sp is the wing planform area. 

The wing is placed on the fuselage according to the 
location of the leading edge of its root chord, determined 
as a fraction of the fuselage length. The distance from the 
nose to the leading edge of the wing is 

x le =1 b r le ( 14 ) 


where Rpp is the ratio of leading edge station to fuselage 
length. 

The first step in computing the wing weight is the 
determination of the geometry of the structural wing box. 
In terms of the input parameters Wjq , ( W/Sp ), aspect ratio 
( AR ), taper ratio (Rtap)> anc * leading edge sweep (A ip), 
the dependent parameters wing area, span, root chord, tip 
chord, and trailing edge wing sweep are computed from 


_ W TO 

w/s 


(15) 


b - ^(AR)S P 


(16) 


C 2-Sp 

* K' + R TAP) 


(17) 


C T = I<TAP C 'r 


(18) 


im(A T E) = tan { A LE) + ~^-{ R TAP- ] ) ( ,9 ) 


(fig. 2). It is assumed that specified portions of the 
streamwise (aerodynamic) chord are required for controls 
and high lift devices, leaving the remainder for the struc- 
tural wing box. The portions of the leading and trailing 
edges that are left for nonstructural use are specified as 
respective fractions Qi and Cy 2 of the streamwise chord. 
Determination of these chord fractions is accomplished 
through visual inspection of the wing planform. Measured 
at the theoretical root chord, the dimensions for the 
leading and trailing edges are 


l LE = C S ] C R 

(20) 

l TE = C S 2 C R 

(21) 


respectively. The intersection of this structural box with 
the fuselage contours determines the location of the 
rectangular carrythrough structure. The width of the 
carrythrough structure, vv£, is defined by the corre- 
sponding fuselage diameter. 

The dimensions of the structural box and of the carry- 
through structure are now determined (fig. 3). The 
structural semispan, is assumed to lie on the quarter- 
chord line, y , whose sweep is given by 

3 1 

tan(A 5 ) = -tan(A L£ ) + -tan(A 7 ' £ ) (22) 


Thus, 

b$ ~ — ~~r~\ (23) 

* 2cos(A 5 ) 

The streamwise chord at any point on the wing is given 
by 

r (Q = c R—^( c R~ c T) (24) 

where £ is measured perpendicular to the vehicle longi- 
tudinal axis from the vehicle centerline toward the 
wingtip. Thus, the streamwise chord is the dimension of 
the wing parallel to the vehicle longitudinal axis. In 
particular, at the wing-fuselage intersection. 


Cr = C' r -j{C' r -C t ) 

(25) 

The structural root and tip chords are 


c SR = ('-c S] -Cs 2 )c r 

(26) 

C ST - (l - c Si ~ C S 2 ) C T 

(27) 
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where Z$(y) and Z(y) are dimensions perpendicular to the 
structural semispan. 

The thickness of the wing box at any spanwise station y is 
determined as a linear interpolation between the root and 
tip thickness ratios multiplied by the chord at y. 

fr/? ; (y), 0 < y < b$ (box structure)] 

'(>’H (32) 

[rR t (0), >’ < 0 (carry through structure)] 

where R((y) is the thickness ratio of the wing as a function 
of position along the quarter chord. 

For the transports in the present study, all the fuel is 
carried within the wing structure. An option is also 
available to carry the fuel entirely within the fuselage, 
negating any bending relief in the wing. (The high 
altitude drone, described in Appendix B, was modeled 
with a fuselage fuel tank.) The volume of the trapezoidal 
planform, rectangular-section wing box structure (includ- 
ing the carrythrough structure) is found as follows: 


respectively. In terms of v, measured along the quarter 
chord from the wing-fuselage intersection toward the 
wingtip, the structural and total chords are given by 

r s{y) = c SR-J-{ c SR~ c ST ) (28) 

b S 

r{y) = C R —^-{C R -C T ) (29) 

b s 

where the structural chord is defined as the dimension of 
the rectangular-section wing box measured parallel to the 
vehicle longitudinal axis. Computation of the widths of 
the wing box and total wing structure, as shown in 
figure 3, is relatively complicated due to the geometry at 
the wingtip and the wing-fuselage intersection. For the 
portion of the wing between the wingtip and the wing- 
fuselage intersection, the respective widths of the wing 
box and total wing structure at any spanwise station yare 

z s{ y) = r s c °s( A s) (30) 

Z(y) = rcos(A 5 ) (31) 
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b S 

V W =2 1 Z s (y)t(y)dy + (l - C S| - C 5z )r [r C 2 rWc 
0 

= 2cos(A5)(l-C 5) -C 52 ) 



0 


x R tR C R -j-[R tR C R -R tT C T ) 
+ { ] ~ C Si ~ c S 2 ) R t R c R w C 

~Q 2 ) cos ( A .g) 


dy 


(33) 


x [R tR C R {2C R + C T ) + R tT C T (C R + 2Cj )] 


+ (*- c 5| ~ c S 2 ) R t R c R w C 

This equation is based on flat upper and lower surfaces 
and neglects the volume taken up by the structure. 


Loads 

Fuselage- Fuselage loading is determined on a station- 
by-station basis along the length of the vehicle. Three 
types of fuselage loads are considered — longitudinal 
acceleration, tank pressure, and bending moment. In the 
present study, all three load types are assumed to occur 
simultaneously to determine maximum compressive and 
tensile loads at the outer shell fibers at each station. 

Bending loads applied to the vehicle fuselage are obtained 
by simulating vehicle pitch-plane motion during a quasi- 
static pull-up maneuver; a landing; and movement over a 
runway bump. Simplified vehicle loading models are 
used where it is assumed that: (1 ) fuselage lift forces 
(nominally zero for subsonic transports) are distributed 
uniformly over the fuselage plan area; (2) wing loading, 
determined independently, is transferred by a couple of 
vertical force and torque through the wing carrythrough 
structure; (3) fuselage weight is distributed uniformly 
over fuselage volume; (4) control surface forces and 
landing gear reactions are point loads; and (5) the pro- 
pulsion system weight, if mounted on the fuselage, is 
uniformly distributed. A factor of safety (nominally 1.5) 
is applied to each load case. The aircraft weight for each 
case is selected as a fraction of gross takeoff weight. The 
resulting one-dimensional loading model is shown in 


figure 4. All fuselage lift forces are assumed to be linear 
functions of angle of attack. Longitudinal bending 
moments are computed for each of the three loading cases 
and the envelope of the maximum values taken as the 
design loading condition. The bending moment computa- 
tion is given in detail in reference 4 and will only be 
summarized here. 

Considering first the pull-up maneuver loading, the 
motion is assumed to be a quasi-static pitch-plane pull-up 
of given normal load factor n (nominally 2.5 for transport 
aircraft). The vehicle is trimmed with the appropriate 
control surface (a horizontal tail for all eight transport 
used for validation in the present study), after which the 
angle of attack is calculated. 

Landing loads are developed as the aircraft descends at a 
given vertical speed, V$, after which it impacts the 
ground; thereafter the main and nose landing gears are 
assumed to exert a constant, or optionally a (1 - cos(coO), 
force during its stroke, S^q, until the aircraft comes to 
rest. The vehicle weight is set equal to the nominal 
landing weight. Wing lift as a fraction of landing weight 
is specified, which reduces the effective load the landing 
gear carries. Likewise, the portion of total vehicle load 
the main gear carries is specified. No pitch-plane motion 
is considered during the landing. 

Runway bump loads are handled by inputting the bump 
load factor into the landing gear. Bump load factor is 
applied according to reference 7. This simulates the 
vehicle running over a bump during taxi. In a similar 
fashion to the landing, the wing lift as a fraction of gross 
takeoff weight is specified, as is the portion of effective 
load input through the main gear. No pitch-plane motion 
is considered during the bump. 

Wing- For the wing, only a quasi-static pull-up maneuver 
condition at load factor n is considered for determining 
loads. At each spanwise station along the quarter chord, 
from the wingtip to the wing-fuselage intersection, the lift 
load, center of pressure, inertia load, center of gravity, 
shear force, and bending moment are computed. For the 
inertia load, it is assumed that the fuel weight Wfj is 
distributed uniformly with respect to the wing volume so 
that the inertial load at y is (W fjW w)*V(y)* where F(y) is 
the volume outboard of y; this volume has centroid Cg(y) 
with respect to station y. An estimate of the wing struc- 
tural weight is included in Wpj for this calculation but the 
calculation is not redone when the actual structural weight 
has been computed. 

There is an option for either a trapezoidal or a Schrenk 
(ref. 8) lift load distribution along the wingspan; the 
trapezoidal distribution represents a uniform lift over the 
wing area (which has a trapezoidal planform) while the 
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Gear Weights 

Figure 4. Loading model. 


Schrenk distribution is an average of the trapezoidal 
distribution with an elliptical distribution, where the lift is 
zero at the wingtip and maximum at the wing-fuselage 
intersection. Prandtl has shown that a true elliptical lift Thus the value of lift at y , Lell* the area of ellipse 

load distribution will have a minimum induced drag, but a outboard of A ell, an d the center of pressure of lift 

combination of the elliptical and trapezoidal distributions outboard of y, Cp EEL , for y measured along the structural 

will give a better representation of actual aircraft loading h° x ma Y determined as 
(ref. 8). 

Plots of trapezoidal and Schrenk lift load distributions are 
shown in figure 5. For the trapezoidal lift load distribution 
the lift load at y is (W/S)Ajpj\p (y), where App^piy) is the 
area outboard of y; the centroid of this area is denoted 

w here y is measured along the quarter chord. 

For the elliptical lift load distribution, the lift load 
matches the contour of an ellipse with the end of its major 
axis on the tip and the end of its minor axis directly above 
the wing-fuselage intersection. The area enclosed by the (36) 

quadrant of the ellipse is set equal to the exposed area of 

the trapezoidal wing panel Cp ( y) = * — (6$ - y) (37) 

LLL 3tc 

respectively. 



(b-wg) 


(l + r tap) c r 
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Elliptical 



For the Schrenk lift load distribution, the average of 
AjRAPiy) and Apppiy) is used to represent the composite 
area, while the average of Cp TRAI Jiy) and Cp EU (y) is 
used to represent the composite center of pressure. 

Using the appropriate outboard area >4 O’) an d center of 
pressure Cp(y), the shear force is 


F s (y) = nK s 


K A -W£L V 
s % 


~- v K 


1=1 


nig 

~2h h ‘g( y igi - y ) w ‘gi 


/= 1 


(38) 


where n e and nig are the number of engines and landing 
gear mounted on the semispan, respectively; W e . and VF/g. 
are the weights of the engine an i^ 1 landing gear, 
respectively; y € . and yig. are the locations of the i { ^ engine 
and i tfl landing gear, respectively; and 


h e[y t -i ->’) = 


i. y e , > y 


o. 


>’e, < >’ 


(39) 




The bending moment is 


yigi > y 

yi 8i <y. 


M{y) = nK s 



Wpj 

Tv 


vc g 


(40) 


n e 

-yj w e i (ye i -y) 

1=1 

n lg 

-XM%, -y) w igi(yi gi ->’) 

i — 1 

(41) 


Structural Analysis 

Fuselage- Weight estimating relationships are now 
developed for the load-carrying fuselage structure. In 
addition, the volume taken up by the fuselage structure is 
also determined. 
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Considering first the circular shell, the stress resultants in 
the axial direction caused by longitudinal bending, axial 
acceleration, and pressure at a fuselage station a are 


N = 
/V XB 

l y 


n *a = 


N *P ~ p 


p 

AP g 


(42) 

(43) 

(44) 


respectively, where r = D/1 is the fuselage radius, A = 7ir 2 
is the fuselage cross-sectional area, and P = 2nr is the 
fuselage perimeter. In equation 42, 1' y = 7tr 3 is the 
moment of inertia of the shell divided by the shell thick- 
ness. In equation 43, for the case of fuselage-mounted 
propulsion, W s is the portion of vehicle weight ahead of 
station jc if a is ahead of the inlet entrance, or the portion 
of vehicle weight behind a if a is behind the nozzle exit. 

In equation 44, Pg is the limit gage pressure differential 
for the passenger compartment during cruise. The total 
tension stress resultant is then 


Nt = N XB +N Xp (45) 

if a is ahead of the nozzle exit, and 

^x=N XB +N Xp+ N XA (46) 

if a is behind it. Similarly, the total compressive stress 
resultant is 


[0, if not pressure stabilized] 


K-N Xb +Nx a 


I N Xp , if stabilized 


if a is ahead of the inlet entrance, and 

( 0, if not pressure stabilized] 

N Xp , if stabilized J 


(47) 


(48) 


if a is behind it. These relations are based on the premise 
that acceleration loads never decrease stress resultants, 
but pressure loads may relieve stress, if pressure stabiliza- 
tion is chosen as an option. The stress resultant in the 
hoop direction is 

N y = rPg K p (49) 

where K p accounts for the fact that not all of the shell 
material (for example, the core material in sandwich 
designs) is available for resisting hoop stress. 


The equivalent isotropic thicknesses of the shell are given 
by 


*S C 



(50) 


h T 


— max(wj, 

Ftu [ 



(51) 


~ Kmg t mg 


(52) 


for designs limited by compressive yield strength (F cv ), 
ultimate tensile strength (F /w ), and minimum gage, 
respectively. In equation 52, t m g is a specified minimum 
material thickness and K mg is a parameter relating to 
t m g which depends on the shell geometry. 

A fourth thickness that must be considered is that for 
buckling critical designs, t $ p , which will now be 
developed. The nominal vehicles of this study have 
integrally stiffened shells stabilized by ring frames. In the 
buckling analysis of these structures, the shell is analyzed 
as a wide column and the frames are sized by the Shanley 
criteria (ref. 6). Expressions are derived for the equivalent 
isotropic thickness of the shell required to preclude 
buckling, tg p , and for the smeared equivalent isotropic 
thickness of the ring frames required to preclude general 
instability, tp . The analysis will be restricted to the case 
of cylindrical shells. The major assumptions are that the 
structural shell behaves as an Euler beam and that all 
structural materials behave elastically. 


For the stiffened shell with frames concept, the common 
procedure of assuming the shell to be a wide column is 
adopted. If the frame spacing is defined as d and Young’s 
modulus of the shell material is defined as £, the buckling 
equation is then 



( h, ^ 


or, solving for t$ B 


% 



(53) 


(54) 


Fuselage structural geometry concepts are presented in 
table 1 ; values of the shell efficiency e for the various 
structural concepts are given in table 2. The structural 
shell geometries available are simply stiffened, 
Z-stiffened, and truss-core sandwich. We next size the 
frames to prevent general instability failure. The Shanley 
criterion is based on the premise that the frames act as 
elastic supports for the wide column; this criterion gives 
the smeared equivalent thickness of the frames as 
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Table 1 . Fuselage structural geometry concepts 


KCON sets 
concept number 

2 Simply stiffened shell, frames, sized for minimum weight in buckling 

3 Z-stiffened shell, frames, best buckling 

4 Z-stiffened shell, frames, buckling-minimum gage compromise 

5 Z-stiffened shell, frames, buckling-pressure compromise 

6 Truss-core sandwich, frames, best buckling 

8 Truss-core sandwich, no frames, best buckling 

9 Truss-core sandwich, no frames, buckling-minimum gage-pressure compromise 


Table 2. Fuselage structural geometry parameters 


Structural concept 
(KCON) 

m 

e 

K mg 

K p 

Kth 

2 

2 

0.656 

2.463 

2.463 

0.0 

3 

2 

0.911 

2.475 

2.475 

0.0 

4 

2 

0.760 

2.039 

1.835 

0.0 

5 

2 

0.760 

2.628 

1.576 

0.0 

6 

2 

0.605 

4.310 

3.965 

0.459 

8 

1.667 

0.4423 

4.820 

3.132 

0.405 

9 

1.667 

0.3615 

3.413 

3.413 

0.320 


tf B 2 r 


nCfN x 

i K F \d 3 E F 


(55) 


where Cp is Shanley’s constant, Kp\ is a frame geometry 
parameter, and Ep is Young’s modulus for the frame 
material. (See ref. 3 for a discussion of the applicability of 
this criterion and for a detailed derivation of the equations 
presented here.) If the structure is buckling critical, the 
total thickness is 


t *S B + l F B 

Minimizing t with respect to d results in 

1 


__ 4 

TiC p 

1 

8 

, 

2 r p F 

K) 

2> 1 

T- 

<N 

1 

{ K F \£ 3 E F E 3 J 

P 

J 


(56) 


(57) 


* S B 



(58) 



d = 


6 r 2 P F I kC pEE 
K . P 1 K F\ E F j 


(59) 


(60) 


where pp is the density of the frame material and p is the 
density of the shell material, so that the shell is three 
times as heavy as the frames. 

Frameless sandwich shell concepts may also be used. 

For these concepts, it is assumed that the elliptical shell 
buckles at the load determined by the maximum compres- 
sive stress resultant N x on the cylinder. The buckling 
equation for these frameless sandwich shell concepts is 
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N 


— = £| 


fz 


rE 


K r J 


( 61 ) 


where m is the buckling equation exponent. Or, solving 
for % 


‘Sr = H 


r N~^ m 


v r£e y 


(62) 


This equation is based on small deflection theory, which 
seems reasonable for sandwich cylindrical shells, 
although it is known to be inaccurate for monocoque 
cylinders. Values of m and e may be found, for example 
in references 9 and 10 for many shell geometries. Table 2 
gives values for sandwich structural concepts available in 
PDCYL, numbers 8 and 9, both of which are truss-core 
sandwich. The quantities N ^ , r, and consequently t $ B , 
will vary with fuselage station dimension x. 

At each fuselage station the shell must satisfy all failure 
criteria and meet all geometric constraints. Thus, the shell 
thickness is selected according to compression, tension, 
minimum gage, and buckling criteria, or 

h = max ( ? % - h T - h G • ‘s B ) (63) 


figure 6, the equivalent isotropic thickness of the smeared 
skin and stringers is 


*$ = '* + 


^fff j b w t w __ 

1 + 1.6 


fivvY 

bs b s 


l b s J 

UJ. 


h 


(67) 


Since only the skin is available for resisting pressure 
loads, 


K 


P 


f 

1 + 1.6 • 

V 


V 


'sA l wJ 


(68) 


For minimum gage designs, if t s > t w then t w = t mg and 


*S = 






(u \ 


+ 1.6 


V b s J 


mg 


(69) 


so that 


K mg ~ 


( . \ 


\‘w j 


+ 1.6 


r b ^ 

u w 

\ b s ) 


(70) 


On the other hand, if t 5 <t w then t s = t mg and 


*S = 


1 + 1.6 

(u \ 
b w 

VT 


b s J 

V f s )_ 


mg 


(71) 


If = t $ B , the structure is buckling critical and the 
equivalent isotropic thickness of the frames, F/r, is 
computed from equation 59. If ig the structure is 

not buckling critical at the optimum frame sizing and the 
frames are resized to make Specifically, a new 

frame spacing is computed from equation 54 as 


d = 


Eei S 2 

Nx 


(64) 


and this value is used in equation 55 to determine tp , 

The total thickness of the fuselage structure is then given 
by the summation of the smeared weights of the shell and 
the frames 


h = h + h (65) 


The shell gage thickness may be computed from 
Fg - Fs / K mg . The ideal fuselage structural weight is 
obtained by summation over the vehicle length 

w ! = 2 (p t s . + pft F . YiAxj (66) 


where the quantities subscripted i depend on x. 

We next discuss the derivation of the structural geometry 
parameters shown in table 2. The Z-stiffened shell, typical 
of modern transport aircraft, will be used as an example 
of skin-stringer-frame construction. Using reference 9 and 


so that 


*m* = « + l-6l 


V ^ 


(72) 


Equations 68, 70, and 72 show that for both pressure 
loading critical and minimum gage limited structure, 
(b w ib s ) and ( t w /t s ) should be as small as possible (he., 
no stringers). As an option in PDCYL, all of the detailed 
shell dimensions shown in figure 6 are computed and 
output at each fuselage station. 

In practice, a typical design will be influenced by bending 
and pressure loads and by the minimum gage constraint, 
and thus a compromise is necessary. If buckling is 
of paramount importance, then a good choice is 
( b w fb s ) = 0.87 and ( t w lt s ) = 1 .06 because this gives the 
maximum buckling efficiency for this concept, namely 
£ = 0.91 1 (ref. 9). From equations 68 and 72, 

K p = K mg = 1 + (1 .6)(0.87)(1 .06) = 2.475 (73) 

This is concept 3 in tables 1 and 2. If pressure dominates 
the loading condition, then ( b w /b s ) =0.6 and (t w /t s )= 0.6 
is a reasonable choice, giving e = 0.76, Kp = 1.576, and 
K mg = 2.628; this is concept 5. For minimum gage 
dominated structure, the geometry (b w lb s ) = 0.58 and 
, (t w lt s ) = 0.90 gives concept 6. 
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Figure 7. Truss-core sandwich geometry. 


The geometry of the truss-core sandwich shell concept e = 0.4423, K m g = 4.820, and Kp = 3. 1 32, concept 8 in 

is shown in figure 7. The equivalent isotropic shell tables 1 and 2. To get a design that is lighter for minimum 

thickness of this concept is gage dominant structure, a geometry is chosen that places 

equal thickness material in the face sheets and the core; 

2 + l _c_ \ t / 74 ) the choice of (tjtf) - 1 .0 and \\f = 45 deg gives structural 

tj cos(\|/) concept 9. These calculations assume that the face sheets 

and core are composed of the same material and are 
Reference 9 shows that the optimum buckling efficiency subject to the same minimum gage constraint, 
is obtained for (tjtf) = 0.65 and \|/ = 55 deg. This gives 
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Since the preceding analysis gives only the ideal weight, 
Wj , the nonoptimum weight, Wyvo (including fasteners, 
cutouts, surface attachments, uniform gage penalties, 
manufacturing constraints, etc.) has yet to be determined. 
The method used will be explained in a later section. 

Wing- Using the geometry and loads applied to the wing 
developed above, the structural dimensions and weight of 
the structural box may now be calculated. The wing struc- 
ture is assumed to be a rectangular multi-web box beam 
with the webs running in the direction of the structural 
semispan. Reference 9 indicates that the critical instability 
mode for multi-web box beams is simultaneous buckling 
of the covers due to local instability and of the webs due 
to flexure induced crushing. This reference gives the 
solidity (ratio of volume of structural material to total 
wing box volume) of the least weight multi-web box 
beams as 


where F5 is the applied shear load and 05 is the allowable 
shear stress. The optimum web spacing (fig. 8) is com- 
puted from (ref. 2) 


%=r 


Mfcj 

(l~ e c)V 2e W 


M 

Zs‘ 2 E; 


2er-3 

) 2 *c 


EC 2e C 


2 eg 


(78) 


where subscripts W and C refer to webs and covers, 
respectively. The equivalent isotropic thicknesses of the 
covers and webs are 


’C~ d W 


M 

Z s tEz c dw 



(79) 


f 

1 = 8 ■ 
V 


M 

Zst 2 E, 


e 


(75) 


where £ and e depend on the cover and web geometries 
(table 3), M is the applied moment, t is the thickness, E is 
the elastic modulus, and Zs is obtained from reference 9. 
The solidity is therefore 


V _ W BENp( y ) 

pz s t 


(76) 


where WpRND IS the weight of bending material per unit 
span and p is the material density. WpRND ls computed 
from equations 75 and 76. The weight per unit span of the 
shear material is 


w shear(>') - 


Pfs 

<*s 


(77) 


t w =t- 
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(80) 


respectively, and the gage thicknesses are 


( 8C ~ K gC* C 

(81) 

’8W ~ K 8W* W 

(82) 


Values of £, e , £<7, Er, and K gc are found in 

table 3 for various structural concepts (ref. 9). If the wing 
structural semispan is divided into N equal length seg- 
ments, the total ideal weight of the wing box structure is 

N 

W BOX = —— T,(W + W SHEARj ) (83) 

1 = 1 


Table 3. Wing structural coefficients and exponents 


Covers 

Webs 

£ 

e 

£ 


£yv 

K 8c 

K 8 W 

Unstiffened 

Truss 

2.25 

0.556 

3.62 

3 

0.605 

1.000 

0.407 

Unstiffened 

Unflanged 

2.21 

0.556 

3.62 

3 

0.656 

1.000 

0.505 

Unstiffened 

Z-stiffened 

2.05 

0.556 

3.62 

3 

0.911 

1.000 

0.405 

Truss 

Truss 

2.44 

0.600 

1.108 

2 

0.605 

0.546 

0.407 

Truss 

Unflanged 

2.40 

0.600 

1.108 

2 

0.656 

0.546 

0.505 

Truss 

Z-stiffened 

2.25 

0.600 

1.108 

2 

0.911 

0.546 

0.405 
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Figure 8. Wing structural concept. 


The wing carrythrough structure consists of torsion 
material in addition to bending and shear material. The 
torsion material is required to resist the twist induced 
due to the sweep of the wing. The bending material is 
computed in a similar manner as that of the box except 
that only the longitudinal component of the bending 
moment contributes. Letting tQ = t(y = 0) and 
M 0 = M(y = 0), 


l c =e 

\ 


Mq cos(As) 
{ 0 C SR e 


V 


) 


(84) 


The weight of the bending material is then 

W BEND C -P^C C SRh) w C < 85 ) 

where wc is the width of the carrythrough structure. 
(When the wing-fuselage intersection occurs entirely 
within the cylindrical midsection, as is the case with all 
eight transport used for validation in the present study, 
w'C = D.) The quantities dw* » ^nd tc are computed in 


the same manner as for the box. The weight of the shear 
material is 


W SHEAR C 



W C 


( 86 ) 


where Fs 0 = F$( 0). 

The torque on the carrythrough structure is 

T=M 0 sin(As) (87) 


and the weight of the torsion material is then 


w TORSION c 


prjtQ -^£srFc 

! o c sr g s 


( 88 ) 


Finally, the ideal weight of the carrythrough structure is 
computed from a summation of the bending shear and 
torsion material, or 

W C = W BEND C + W SHEAR C + w TORSION c ( 89 ) 
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As in the case of the fuselage structural weight, nonopti- 
mum weight must be added to the ideal weight to obtain 
the true wing structural weight. The method used will be 
discussed below. 

The static deflection of the wingtip under the pull-up 
maneuver is also determined. Using the moment-area 
method applied to an Euler beam (ref. 1 1 ), the deviation 
of point B on the deflected surface from the tangent 
drawn from another point A on the surface is equal to the 
area under the M/(Ef) diagram between A and B multi- 
plied by the distance to the centroid of this area from B , 

A A 

'ba = J* yd& = J ~y d y ( 90 ) 

B B 

where 0 is the angular displacement of the beam and y is 
the longitudinal axis of the beam. For the case of a wing 
with trapezoidal planform, the longitudinal axis, y, will lie 
along the quarter-chord line (fig. 3). For a wing with a 
horizontal unloaded configuration, the tangential devia- 
tion, t£A, will equal the true vertical tip displacement 
(assumed to be the case). Only the wing cover contributes 
to the bending resistance, while the webs offer similar 
shear stiffness. The wing area moment of inertia, /, at 
any structural semispan station y is determined with the 
Parallel Axis theorem, as cover thickness is small when 
compared with total wing thickness. 


Regression Analysis 

Overview— Using fuselage and wing weight statements 
of eight subsonic transports, a relation between the calcu- 
lated load-bearing structure weights obtained through 
PDCYL and the actual load-bearing structure weights, 
primary structure weights, and total weights is determined 
using statistical analysis techniques. A basic application 
which is first described is linear regression, wherein the 
estimated weights of the aircraft are related to the weights 
calculated by PDCYL with a straight line, v = mx + b , 
where y is the value of the estimated weight, m is the 
slope of the line, x is the value obtained through PDCYL, 
and b is the y-intercept. This line is termed a regression 
line, and is found by using the method of least squares, 
in which the sum of the squares of the residual errors 
between actual data points and the corresponding points 
on the regression line is minimized. Effectively, a straight 
line is drawn through a set of ordered pairs of data (in 
this case eight weights obtained through PDCYL and the 
corresponding actual weights) so that the aggregate 
deviation of the actual weights above or below this line is 
minimized. The estimated weight is therefore dependent 
upon the independent PDCYL weight. 


As an example, if the form of the regression equation is 
linear, the estimated weight is 

y est =mx caic +b (91) 

where m is the slope, b is the intercept, and x ca [ c is the 
weight PDCYL calculates. The resulting residual to be 
minimized is 


n 2 

^ ^ {y actual i ~ y'estj ) (92) 

i=l 


or 

n 2 

^ actual i “ nix calcj ~ (93) 

/=1 

where y actual is the actual component weight and/? is the 
number of aircraft whose data are to be used in the fit. By 
taking partial derivatives of the residual error with respect 
to both m and b, equations for the values of these two 
unknown variables are found to be 


ii ii n 

^ Xca ^ c i y act j ~ x calc i y act , 


m — - 


i~\ 


i=l 


j= 1 


ll u 

"Tx-ala ~ Y x calcj 


i= 1 

h = y ( icl ~ "Xalc 


W= 1 

.v, v = mean values of a* and v 


(94) 


(95) 


Of key importance is the degree of accuracy to which the 
prediction techniques are able to estimate actual aircraft 
weight. A measure of this accuracy, the correlation 
coefficient, denoted R, represents the reduction in residual 
error due to the regression technique. R is defined as 


R : 


E r 


(96) 


where E t and E r refer to the residual errors associated 
with the regression before and after analysis is performed, 
respectively. A value of R = I denotes a perfect fit of the 
data with the regression line. Conversely, a value of R = 0 
denotes no improvement in the data fit due to regression 
analysis. 

There are two basic forms of equations which are 
implemented in this study. The first is of the form 

y'est ~ mx calc (97) 

The second general form is 


y'est 


= »<„c 


(98) 
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The first form is a simplified version of the linear 
example as discussed above, with the v-intercept term set 
to zero. However, because the second general equation is 
not linear, nor can it be transformed to a linear equation, 
an alternative method must be employed. In order to 
formulate the resulting power-intercept regression equa- 
tion, an iterative approach developed by D. W. Marquardt 
is utilized (ref. 12). This algorithm starts at a certain point 
in space, and, by applying the method of steepest descent, 
a gradient is obtained which indicates the direction in 
which the most rapid decrease in the residual errors will 
occur. In addition, the Taylor Series method produces a 
second similar vector. Interpolation between these two 
vectors yields a direction in which to move the point in 
order to minimize the associated error. After several 
iterations, the process converges to a minimum value. It 
should be noted that there may be several local minimums 
and there is no guarantee that the method converges to the 
global one. 

Fuselage- The analysis above is used to develop a 
relationship between weight calculated by PDCYL and 
actual wing and fuselage weights. The data were obtained 
from detailed weight breakdowns of eight transport air- 
craft (refs. 13-17) and are shown in table 4 for the 
fuselage. Because the theory used in the PDCYL analysis 
only predicts the load-carrying structure of the aircraft 
components, a correlation between the predicted weight 
and the actual load-carrying structural weight and primary 
weight, as well as the total weight of the fuselage, was 
made. 


Structural weight consists of all load-carrying members 
including bulkheads and frames, minor frames, covering, 
covering stiffeners, and longerons. For the linear curve- 
fit, the resulting regression equation is 

^actual = 1 3503 W calc R = 0.9946 (99) 

This shows that the nonoptimum factor for fuselage 
structure is 1 .3503; in other words, the calculated weight 
must be increased by about 35 percent to get the actual 
structural weight. For the alternative power- intercept 
curve fitting analysis, the resulting load-carrying 
regression equation is 

Kauai = 1.1 304^ 79 R = 0.9946 (100) 

To use either of these equations to estimate total fuselage 
weight, nonstructural weight items must be estimated 
independently and added to the structural weight. 

Primary weight consists of all load-carrying members as 
well as any secondary structural items such as joints 
fasteners, keel beam, fail-safe straps, flooring, flooring 
structural supplies, and pressure web. It also includes the 
lavatory structure, galley support, partitions, shear ties, tie 
rods, structural firewall, torque boxes, and attachment 
fittings. The linear curve fit for this weight yields the 
following primary regression equation 

Kauai = 1 -8872 W calc R = 0.99 17 (101) 

The primary power- intercept regression equation is 

Kauai = 1 6399 W caU- 4 ’ K = 0.9917 (102) 


Table 4. Fuselage weight breakdowns for eight transport aircraft 


Weight, lb 

Aircraft 

PDCYL 

Load-carrying structure Primary structure 

Total structure 

B-720 

6545 

9013 

13336 

19383 

B-727 

5888 

8790 

12424 

17586 

B-737 

3428 

5089 

7435 

11831 

B-747 

28039 

39936 

55207 

72659 

DC-8 

9527 

13312 

18584 

24886 

MD-1 1 

20915 

25970 

34999 

54936 

MD-83 

7443 

9410 

11880 

16432 

L- 1011 

21608 

28352 

41804 

52329 
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The total fuselage weight accounts for all members of the 
body, including the structural weight and primary weight. 
It does not include passenger accommodations, such as 
seats, lavatories, kitchens, stowage, and lighting; the 
electrical system; flight and navigation systems; alighting 
gear; fuel and propulsion systems; hydraulic and pneu- 
matic systems; the communication system; cargo accom- 
modations; flight deck accommodations; air conditioning 
equipment; the auxiliary power system; and emergency 
systems. Linear regression results in the following total 
fuselage weight equation 

Wactual = 2.5686 W cdc R = 0.9944 (103) 

This shows that the nonoptimum factor for the total 
fuselage weight is 2.5686; in other words, the fuselage 


structure weight estimated by PDCYL must be increased 
by about 157 percent to get the actual total fuselage 
weight This nonoptimum factor is used to compare 
fuselage structure weight estimates from PDCYL with 
total fuselage weight estimates from the Sanders and the 
Air Force equations used by ACSYNT. 

The total fuselage weight power-intercept regression 
equation is 

W actual = 3.9089W C ^ 578 R = 0.9949 (104) 

Plots of actual fuselage component weight versus 
PDCYL-calculated weight, as well as the corresponding 
linear regressions, are shown in figures 9-11. 
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Figure 9. Fuselage load-carrying structure and linear regression. 
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Table 5. Wing weight breakdowns for eight transport aircraft 


Weight, lb 

Aircraft 

PDCYL 

Load-carrying structure 

Primary structure 

Total structure 

B-720 

13962 

11747 


18914 

23528 

B-727 

8688 

8791 


12388 

17860 

B-737 

5717 

5414 


7671 

10687 

B-747 

52950 

50395 


68761 

88202 

DC-8 

22080 

19130 


27924 

35330 

MD-1 1 

33617 

35157 


47614 

62985 

MD-83 

6953 

8720 


1 1553 

15839 

L-101 1 

25034 

28355 


36101 

46233 


Wing- The same analysis was performed on the wing 
weight for the sample aircraft and is shown in table 5. 

The wing box, or load-carrying structure, consists of spar 
caps, interspar coverings, spanwise stiffeners, spar webs, 
spar stiffeners, and interspar ribs. The wing box linear 
regression equation is 

Wactuai =0.9843 W cak /? = 0.9898 (105) 

so that the nonoptimum factor is 0.9843. Power-intercept 
regression results in 

w acllla , = 1 .3342 1 R = 0.9902 (106) 

Wing primary structural weight includes all wing box 
items in addition to auxiliary spar caps and spar webs, 
joints and fasteners, landing gear support beam, leading 
and trailing edges, tips, structural firewall, bulkheads, 
jacket fittings, terminal fittings, and attachments. Linear 
regression results in 

Kauai = 1 -3442 W calc R = 0.9958 (107) 

Power-intercept regression yields 

w actual = 2 . 1 926 W°& 534 R = 0.9969 (108) 


The total wing weight includes wing box and primary 
weight items in addition to high-lift devices, control 
surfaces, and access items. It does not include the pro- 
pulsion system, fuel system, and thrust reversers; the 
electrical system; alighting gear; hydraulic and pneumatic 
systems; anti-icing devices; and emergency systems. The 
resulting total weight linear regression equation is 

^actual = 1 -7372 W calc R = 0.9925 ( 1 09) 

This shows that the nonoptimum factor for the total wing 
weight is 1 .7372; in other words, the wing box weight 
estimated by PDCYL must be increased by about 
74 percent to get the actual total wing weight. This 
nonoptimum factor is used to compare wing box weight 
estimates from PDCYL with total wing weight estimates 
from the Sanders and the Air Force equations used by 
ACSYNT. 

The power-intercept equation for total wing weight is 

Wactual = 3-7464 W°;£ 268 R = 0.9946 (110) 

Plots of actual wing component weight versus PDCYL- 
calculated weight, as well as the corresponding linear 
regressions, are shown in figures 12-14. 
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Figure 14. Wing total structure and linear regression. 


Discussion- Both fuselage and wing weight linear and 
power regressions give excellent correlation with the 
respective weights of existing aircraft, as evidenced by 
the high values of the correlation coefficient, R. It should 
be noted that even though the power-based regressions 
give correlations equal to or better than the linear regres- 
sions their factors may vary distinctly from the linear 
cases. This is due to their powers not equaling unity. 


Because estimates of non-load-bearing primary structure 
are generally not available at the conceptual design stage, 
and because nonprimary structure is probably not well 
estimated by a nonoptimum factor, equations 101 and 107 
are recommended for estimating the primary structural 
weights of the respective transport fuselage and wing 
structures (figs. 10 and 13). 


24 




Appendix A - User’s Manual, Example 

Description 

The purpose of this appendix is to give a detailed example 
of the input procedure used to allow PDCYL to calculate 
fuselage and wing weights for a sample transport aircraft 
during an ACSYNT run. A sample output from PDCYL 
will also be given. The Boeing 747-2 IP will be used for 
the example. The layout of the 747-2 IP is shown in 
figure 15. The weights of the load-carrying portions of the 
fuselage and wing box for the 747-2 IP will be calculated 
by PDCYL and scaled by the respective nonoptimum 
factors developed earlier to give estimates for the weights 
of the fuselage and wing. A comparison between methods 
currently used by ACSYNT to estimate fuselage and wing 
weights and PDCYL output will be made with the 
corresponding actual weights of the 747-2 IP. 

Input 

PDCYL requires input from both the existing ACSYNT 
data structure and an additional namelist containing data 
required by PDCYL which are not contained within the 
current ACSYNT format. There are three steps to run 
PDCYL within ACSYNT. First, the aircraft type is 
specified in the ACSYNT Control input. Currently the 
Transport Aircraft type is used. Second, data within 
ACSYNT module namelists are required. The ACSYNT 
Geometry, Trajectory, and Weights modules supply data 
for PDCYL execution. PDCYL uses the WING, HTAIL, 
VTAIL, FUS, WPOD, and FPOD namelists from the 
Geometry module. From the Trajectory module, the 
TRDATA namelist is used. From the Weights module the 
OPTS namelist is used. Third, data from the PDCYLIN 
namelist are used. 

Variables used from ACSYNT namelists and the 
PDCYLIN namelist are given in tables 6 and 7, respec- 
tively. Default values for all variables are also given. 
These default values match the Boeing 747-2 IP. Key 
configuration parameters are given for each of the eight 
aircraft used in the validation study in table 8. An 
example of the PDCYLIN namelist input for the 747-2 IP 
is shown in figure 16. 

A description of the specific structural concepts used to 
model both the fuselage and wing is given in the Struc- 
tural Analysis section. As was noted earlier, the typical 
modern transport aircraft fuselage is a Z-stiffened shell. 
The buckling-minimum material gage compromise was 


employed because it gives the Iowest-weight (optimal) 
structure for the eight aircraft investigated in this study. 

Output 

PDCYL weights output begins with the wing box and 
carrythrough structure analysis. The wing is sized during 
a quasi-static pull-up maneuver where the load factor is 
set equal to the ultimate load factor (nominally 3.75). 
Wing output contains three parts. First is the overall 
geometrical configuration. Second is a detailed station- 
by-station bending, shear, and torsion analysis and 
corresponding geometrical sizing along the span. Third 
is the detailed geometrical layout, loading, and weight 
breakdown of the carrythrough structure, weight break- 
down of the wing components, and deflection of the 
wingtip. This wing weight is multiplied by the nonopti- 
mum factor and returned to ACSYNT. An example of the 
PDCYL wing weight output for the 747-2 IP is shown in 
figure 17. 

Next, the fuselage is analyzed. Fuselage output contains 
four parts. First is the overall geometrical layout and 
weight breakdown. Second is a station-by-station bend- 
ing, shear, and axial stress analysis. Up to three load cases 
are investigated. In order they are a quasi-static pull-up 
maneuver, a landing maneuver, and travel over runway 
bumps. Third, the envelope of worst-case loading is 
shown for each station, from which the shell and frames 
are sized. Corresponding unit weight breakdowns are also 
given. As an option, the detailed geometric configuration 
at each station may be output. Fourth, weights summaries 
are given for the top and bottom sections of the fuselage 
(nominally the same). These summaries are then averaged 
to give the weight summary of the entire fuselage" The 
fuselage weight, including the corresponding nonopti- 
mum factor, is returned to ACSYNT. An example of the 
PDCYL fuselage weight output for the 747-2 IP is shown 
in figure 18. 

Figure 19(a) shows a comparison between fuselage 
weight estimates from the Sanders equation, the Air Force 
equation, and PDCYL with the actual fuselage weight of 
the 747-2 IP. Figure 19(b) shows a similar comparison for 
the wing weight. SLOPE and TECH factors were set to 
one for the comparisons in Figures 19(a) and 19(b), while 
the nonoptimum factors are those relating PDCYL esti- 
mations of structure weight to respective total component 
weight. 
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Table 6. ACSYNT variables 


Variable 

Type 

Dimension 

Description 

Units/comment 

Default (747) 

1. Geometry module 

Namelist WING 

SWEEP 

float 

1 

Sweep of wing. 

degrees 

37.17 

KSWEEP 

integer 

1 

1 —» Referenced to the leading edge. 

2 — » Referenced to the quarter chord. 

3 — > Referenced to the trailing edge. 


2 

AR 

float 

1 

Aspect ratio of wing. 


6.96 

TAPER 

float 

1 

Taper ratio of wing. 


0.2646 

TCROOT 

float 

1 

Thickness-to-chord ratio at the root. 


0.1794 

TCTIP 

float 

1 

Thickness-to-chord ratio at the tip. 


0.078 

ZROOT 

float 

1 

Elevation of MAC above fuselage 


-0.1 




reference plane, measured as a fraction 
of the local fuselage radius. 



AREA 

float 

1 

Planform area of wing. 

ft 2 

5469 

DIHED 

float 

1 

Dihedral angle of wing. 

degrees 

7 

XWING 

float 

1 

Ratio of distance measured from nose to 


0.249 




leading edge of wing to total fuselage 
length. 



Namelist HTAIL (horizontal tail) 

SWEEP 

float 

1 

Sweep of tail 

degrees 

34.29 

KSWEEP 

integer 

1 

1 — » Referenced to the leading edge. 

2 — > Referenced to the quarter chord. 

3 -> Referenced to the trailing edge. 


2 

AR 

float 

1 

Aspect ratio of the horizontal wing. 

(span) 2 /area 

3.625 

TAPER 

float 

1 

Taper ratio of the horizontal wing. 

tip chord/root 

0.25 





chord 


TCROOT 

float 

1 

Thickness-to-chord ratio at the root. 


0.11 

TCTIP 

float 

1 

Thickness-to-chord ratio at the tip. 


0.08 

ZROOT 

float 

1 

Elevation of MAC above fuselage 


0.69 




reference plane, measured as a fraction 
of the local fuselage radius. 



AREA 

float 

1 

Planform area of the horizontal wing. 

ft 2 

1470 

XHTAIL 

float 

1 

Position for trailing edge of tail root 


1 




chord. If ZROOT < 1, then XHTAIL is 
given as a fraction of body length. Else, 
XHTAIL is given as a fraction of the 
local vertical tail chord. 




27 


Table 6. Continued 


Variable 

Type Dimension 

Description 

Units/comment 

Default (747) 

Namelist VTAIL (vertical tail) 

SWEEP 

float 1 

Sweep of vertical tail. 

degrees 

45.73 

KSWEEP 

integer 1 

1 — » Referenced to the leading edge. 

2 — » Referenced to the quarter chord. 

3 — » Referenced t the trailing edge. 


2 

AR 

float 1 

Aspect ratio of vertical tail. 

(span) 2 /area 

1.247 

TAPER 

float 1 

Taper ratio of vertical tail. 

tip chord/root 

0.34 




chord 


TCROOT 

float 1 

Thickness-to-chord ratio at root. 


0.1298 

TCTIP 

float 1 

Thickness-to-chord ratio at tip. 


0.089 

ZROOT 

float 1 

Elevation of MAC above fuselage 


0.6 



reference plane, measured as a fraction 
of the local fuselage radius. 



AREA 

float 1 

Planform area of vertical tail. 

it* 

830 

Namelist FUS (fuselage) 

FRN 

float 1 

Fineness ratio of the nose section. 

length/diameter 

2.13 

FRAB 

float 1 

Fineness ratio of after-body section. 

length/diameter 

3.29 

BODL 

float 1 

Length of fuselage. 

ft 

225.167 

BDMAX 

float 1 

Maximum diameter of fuselage. 

ft 

20.2 


Namelist WPOD (wing-mounted propulsion pod) 


DIAM 

float 

1 Engine diameter. 

ft 

6.2 

LENGTH 

float 

1 Length of engine pod. 

ft 

15 

X 

float 

1 X location of nose of pod relative to 

leading edge of wing, given as a 
fraction of local chord of wing (>0 if 
face of pod is behind leading edge of 
wing). 


-0.631 

Y 

float 

1 Y location of center of pod, given as a 

fraction of semispan, measured from 
body centerline. 


0.241 

Z 

float 

1 Z location of center of pod above wing 

local chord, given as fraction of 
maximum pod diameter. 


-0.83 

SWFACT 

float 

1 Wetted area multiplier. 
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Table 6. Concluded 


Variable 

Type 

Dimension 

Description 

Units/comment 

Default (747) 

Namelist FPOD (fuselage-mounted propulsion pod) 

DIAM 

float 

1 

Engine diameter. 

ft 

N/A 

LENGTH 

float 

i 

Length of engine pod. 

ft 

N/A 

SOD 

float 

1 

Stand-off-distance, the distance from the 
pod wall to the fuselage wall, given as a 
fraction of maximum pod radius. 


N/A 

THETA 

float 

1 

Angular orientation of pod, THETA 
measured positive up from the 
horizontal reference plane. 

degrees 

N/A 

X 

float 

1 

X location of nose relative to nose of 
fuselage, given as a fraction of body 
length. 


N/A 


2. Trajectory module 


Namelist TRDATA (used for load factors) 


DESLF float 

ULTLF float 

1 Design load factor. 

1 Ultimate load factor, usually 1 .5*DESLF. 

N/A 

N/A 

2.5 

3.75 

3. Weights module 




Namelist OPTS 

WGTO float 

1 Gross take-off weight. 

lb 

713000 

WE float 

1 Total weight of propulsion system 

(includes both wing and fuselage 
mounted engines). 

lb 

44290 
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Table 7. PDCYL variables 


Variable 

Type 

Dimension 

Description 

Units/commeni 

Default (747) 

Namelist PDCYLIN 





Wing 






Material properties 

PS 

float 

1 

Plasticity factor. 


1 

TMGW 

float 

1 

Min. gage thickness for the wing 

inches 

0.2 

EFFW 

float 

1 

Buckling efficiency of the web. 


0.656 

EFFC 

float 

1 

Buckling efficiency of the covers. 


1.03 

ESW 

float 

1 

Young’s Modulus for wing material. 

psi 

1.07E+07 

FCSW 

float 

1 

Ult. compressive strength of wing. 

psi 

54000 

DSW 

float 

1 

Density of the wing material. 

lb/in. 3 

0.101 

KDEW 

float 

1 

Knock-down factor for Young's Modulus. 


1 

KDFW 

float 

1 

Knock-down factor for Ultimate strength. 


1 

Geometric parameters 

ISTAMA 

integer 

1 

1 the position of the wing is unknown. 

2 the position of the wing is known. 


2 

CS1 

float 

1 

Position of structural wing box from 
leading edge as percent of root chord. 


0.088 

CS2 

float 

1 

Position of structural wing box from 
trailing edge as percent of root chord. 


0.277 

Structural concept 

CLAQR 

float 

1 

Ratio of body lift to wing lift. 

For subsonic 
aircraft 
CLAQR ~ 0.0 

0.001 

I FUEL 

integer 

1 

1 — > no fuel is stored in the wing. 

2 — > fuel is stored in the wing. 


2 

CWMAN 

float 

1 

Design maneuver load factor. 


1 

CF 



Shanley's const, for frame bending. 


6.25E-05 


Fuselage 


Structural concept 


CKF 

float 

1 

Frame stiffness coefficient. 

5.24 

EC 

float 

1 

Power in approximation equation for 
buckling stability. 

2.36 

KGC 

float 

1 

Buckling coefficient for component 
general buckling of stiffener web panel. 

0.368 

KGW 

float 

1 

Buckling coefficient for component local 
buckling of web panel. 

0.505 
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Table 7. Continued 


KCON(T/B) Structural geometry concept Default (747) 

2 Simply stiffened shell, frames, sized for minimum weight in buckling 

3 Z-stiffened shell, frames, best buckling 

4 Z-stiffened shell, frames, buckling-minimum gage compromise 

5 Z-stiffened shell, frames, buckling-pressure compromise 4 

6 Truss-core sandwich, frames, best buckling 

8 Truss-core sandwich, no frames, best buckling 

9 Truss-core sandwich, no frames, buckling-min. gage-pressure compromise 


Variable 

Type 

Dimension 

Description 

Units/Comment 

Default (747) 

Material properties 

FTS(T/B) 

float 

4 

Tensile strength on (top/bottom). 

psi 

58500 

FCS(T/B) 

float 

4 

Compressive strength. 

psi 

54000 

ES(T/B) 

float 

4 

Young’s Modulus for the shells. 

psi 

1.07E+07 

EF(T/B) 

float 

4 

Young's Modulus for the frames. 

psi 

1 .07E+07 

DS(T/B) 

float 

4 

Density of shell material on (t/b). 

Ib/in. 3 

0.101 

DF(T/B) 

float 

4 

Density of frame material. 

Ib/in. 3 

0.101 

TMG(T/B) 

float 

4 

Minimum gage thickness. 

in. 

0.071 

KDE 

float 

1 

Knock-down factor for modulus. 


1 

KDF 

float 

1 

Knock-down factor for strength. 


1 

Geometric parameters 

CLBR1 

float 

1 

Fuselage break point as a fraction of total 


1.1 




fuselage length. 



ICYL 

integer 

1 

1 — > modeled with a mid-body cylinder. 
Else — > use two power-law bodies back to 


1 




back. 



Loads 

AXAC 

float 

1 

Axial acceleration. 

g’s 

0 

CMAN 

float 

1 

Weight fraction at maneuver. 


1 

ILOAD 

integer 

1 

1 — » analyze maneuver only. 

2 — > analyze maneuver and landing only. 

3 — > analyze bump, landing, and 


3 




maneuver. 



PG(T/B) 

float 

12 

Fuselage gage pressure on (top/bot). 

psi 

13.65 

WFBUMP 

float 

1 

Weight fraction at bump. 


0.001 

WFLAND 

float 

1 

Weight fraction at landing. 


0.9 
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Table 7. Continued 


Variable 

Type 

Dimension 

Description 

Units/Comment 

Default (747) 

Landing gear 

VSINK 

float 

i 

Design sink velocity at landing. 

ft/sec 

10 

STROKE 

float 

1 

Stroke of landing gear. 

ft 

2.21 

CLRG1 

float 

1 

Length fraction of nose landing gear. 


0.1131 

CLRG2 

float 

1 

Length fraction of main landing gear 


0.466 




measured as a fraction of total fuselage 
length. 



WFGR1 

float 

1 

Weight fraction of nose landing gear. 


0.0047 

WFGR2 

float 

1 

Weight fraction of main landing gear. 


0.0398 

IGEAR 

integer 

1 

1 — » main landing gear located on 
fuselage. 

2 main landing gear located on wing. 


2 

GFRL 

float 

1 

Ratio of force taken by nose landing gear 


0.001 




to force taken by main gear at landing. 



CLRGW1 

float 

1 

Position of wing gear as a fraction of 

If only 1 wing 

0.064 

CLRGW2 

float 

1 

structural semispan. 

gear, set 

0.1844 





CLRGW2 = 0.0 


Tails 

IT AIL 

integer 

1 

1 — » control surfaces mounted on tail. 


1 




2 — > control surfaces mounted on wing. 



Weights 

WTFF 

float 

1 

Weight fraction of fuel. 


0.262 

CBUM 

float 

1 

Weight fraction at bump. 


1 

CLAN 

float 

1 

Weight fraction at landing. 


0.791 
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Table 7. Concluded 


Variable 

Type 

Dimension 

Description Units/Comment 

Default (747) 

Factors 





ISCHRENK 

integer 

1 

1 -4 use Schrenk load distribution on 
wing. 

Else -4 use trapezoidal distribution. 

1 

ICOMND 

integer 

1 

1 -4 print gross shell dimensions 
envelope. 

2 -4 print detailed shell geometry. 

1 

WGNO 

float 

1 

Nonoptimal factor for wing (including the 
secondary structure). 

i 

SLFMB 

float 

1 

Static load factor for bumps. 

1.2 

WMIS 

float 

1 

Volume component of secondary 
structure. 

0 

WSUR 

float 

1 

Surface area component of secondary 
structure. 

0 

WCW 

float 

1 

Factor in weight equation for nonoptimal 
weights. 

1 

WCA 

float 

1 

Factor in weight equation multiplying 
surface areas for nonoptimal weights. 

0 

NWING 

integer 

1 

Number of wing segments for analysis. 
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Table 8. Key configuration parameters for eight transport aircraft 


Variable 720 727 737 747 DC-8 MD-11 MD-83 L-1011 


ACSYNT INPUT PARAMETERS 


1. Geometry module 


Namelist WING 


SWEEP 

35 

32 

25 

37.17 

30.6 

35 

24.16 

35 

KSWEEP 

2 

2 

2 

2 

2 

2 

2 

2 

AR 

6.958 

7.67 

8.21 

6.96 

7.52 

7.5 

9.62 

6.98 

TAPER 

0.333 

0.2646 

0.2197 

0.2646 

0.1974 

0.255 

0.156 

0.3 

TCROOT 

0.1551 

0.154 

0.126 

0.1794 

0.1256 

0.167 

0.138 

0.13 

TCTIP 

0.0902 

0.09 

0.112 

0.078 

0.105 

0.093 

0.12 

0.09 

ZROOT 

-1 

-1 

-0.25 

-0.1 

-1 

-0.79 

-1 

-1 

AREA 

2460 

1587 

1005 

5469 

2927 

3648 

1270 

3590 

DIHED 

3 

3 

6 

7 

3 

6 

3 

3 

XWING 

0.2963 

0.376 

0.35 

0.249 

0.302 

0.218 

0.468 

0.359 


Namelist HTAIL 


SWEEP 

35 

31.05 

30.298 

34.29 

35 

35.5 

30.8 

3.5 

KSWEEP 

2 

2 

2 

2 

2 

2 

2 

2 

AR 

3.15 

3.4 

4.04 

3.625 

4.04 

3.43 

4.88 

4 

TAPER 

0.457 

0.383 

0.3974 

0.25 

0.329 

0.412 

0.357 

0.33 

TCROOT 

0.1 1 

0.11 

0.132 

0.11 

0.095 

0.143 

0.107 

0.095 

TCTIP 

0.09 

0.0894 

0,108 

0.08 

0.08 

0.1067 

0.08 

0.08 

ZROOT 

0.5 

2 

0.67 

0.69 

0.25 

0.6875 

2 

0.5 

AREA 

500 

376 

312 

1470 

559 

920 

314 

1282 

XHTAIL 

1 

0.95 

0.8532 

0.974 

1 

0.96 

0.98 

0.9265 


Namelist VTAIL 


SWEEP 

35 

48.4 

34.16 

45.73 

35 

38 

39.4 

35 

KSWEEP 

2 

2 

2 

2 

2 

2 

2 

2 

AR 

1.45 

1.09 

1.814 

1.247 

1.905 

1.73 

1.48 

1.6 

TAPER 

0.484 

0.641 

0.3024 

0.34 

0.292 

0.343 

0.844 

0.3 

TCROOT 

0.1 1 

0.1 1 

0.1322 

0.1298 

0.096 

0.105 

0.127 

0.1 1 

TCTIP 

0.0896 

0.09 

0.1081 

0.089 

0.101 

0.125 

0.103 

0.0896 

ZROOT 

0.95 

0.2 

0 

0.6 

0.95 

0.85 

0.9 

0.95 

AREA 

312.4 

356 

225 

830 

352 

605 

550 

550 
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Table 8. Continued 


Variable 720 727 737 747 DC-8 MD-11 MD-83 L-1011 


Namelist FUS 

FRN 

1.81 

2 

1.915 

2.13 

2 

1.67 

1.15 

1.76 

FRAB 

2.86 

2.831 

2.361 

3.29 

2.9375 

2.27 

2.73 

2.96 

BODL 

130.5 

116.67 

90.58 

225.167 

153 

192.42 

135.5 

177.67 

BDMAX 

14.21 

14.2 

13.167 

20.2 

13.5 

19.75 

1 1.44 

19.583 

Namelist WPOD (inboard) 

DIAM 

3.24 

N/A 

3.542 

6.2 

4.42 

9.04 

N/A 

3.24 

LENGTH 

12.15 

N/A 

10 

15 

12.15 

18.08 

N/A 

12.15 

X 

0.917 

N/A 

-0.22 

-0.631 

-0.4 

-0.558 

N/A 

-0.639 

Y 

0.386 

N/A 

0.343 

0.241 

0.352 

0.33125 

N/A 

0.461 

Z 

-1 

N/A 

-0.548 

-0.83 

-1.2 

-0.5 

N/A 

-1 

SWFACT 

1 

N/A 

1 

1 

1 

l 

N/A 

1 

Namelist WPOD (outboard) 

DIAM 

3.24 

N/A 

N/A 

6.2 

4.42 

N/A 

N/A 

N/A 

LENGTH 

12.15 

N/A 

N/A 

15 

12.15 

N/A 

N/A 

N/A 

X 

0.917 

N/A 

N/A 

-0.631 

-0.955 

N/A 

N/A 

N/A 

Y 

0.674 

N/A 

N/A 

0.441 

0.61 

N/A 

N/A 

N/A 

Z 

-1 

N/A 

N/A 

-0.83 

-1.2 

N/A 

N/A 

N/A 

SWFACT 

1 

N/A 

N/A 

1 

1 

N/A 

N/A 

N/A 

Namelist FPOD 

DIAM 

N/A 

3.542 

N/A 

N/A 

N/A 

9.04 

6.6 

3.24 

LENGTH 

N/A 

10 

N/A 

N/A 

N/A 

40.68 

20.34 

12.15 

SOD 

N/A 

0 

N/A 

N/A 

N/A 

0 

0 

0 

THETA 

N/A 

90 

N/A 

N/A 

N/A 

90 

0 

90 

X 

N/A 

0.699 

N/A 

N/A 

N/A 

0.812 

0.746 

0.725 

SYMCOD 

N/A 

1 

N/A 

N/A 

N/A 

1 

0 

-1 

Namelist FPOD (third engine) 

DIAM 

N/A 

3.542 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

LENGTH 

N/A 

10 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

SOD 

N/A 

0.2 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

THETA 

N/A 

14.8 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

X 

N/A 

0.699 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 

SYMCOD 

N/A 

0 

N/A 

N/A 

N/A 

N/A 

N/A 

N/A 


Table 8. Continued 


Variable 720 

727 

737 

747 

DC-8 

MD-11 

MD-83 

L-1011 

2. Trajectory module 








Namelist TRDATA 

DESLF 2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

ULTLF 3.75 

3.75 

3.75 

3.75 

3.75 

3.75 

3.75 

3.75 

3. Weights module 








Namelist OPTS 

WGTO 202000 

160000 

100800 

713000 

335000 

602500 

140000 

409000 

Namelist FIXW 








WE 18202 

12759 

8165 

44290 

27058 

40955 

10340 

34797 


PDCYL INPUT PARAMETERS 


Wing 


Geometric parameters 

ISTAMA 

2 

2 

2 

2 

2 

2 

2 

2 

CS1 

0.1 

0.2125 

0.0724 

0.088 

0.0818 

0.168 

0.181 

0.093 

CS2 

0.27 

0.25 

0.238 

0.277 

0.136 

0.2835 

0.271 

0.296 

Structural concept 

CLAQR 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

IFUEL 

2 

2 

2 

2 

2 

2 

2 

2 

CWMAN 

1 

1 

I 

1 

1 

I 

1 

1 

CF 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

6.25E-05 

Material properties 

PS 

1 

1 

1 

1 

I 

1 

1 

1 

TMGW 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

EFFW 

0.656 

0.656 

0.656 

0.656 

0.656 

0.656 

0.656 

0.656 

EFFC 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

1.03 

ESW 

1.08E+07 

1 .08E+07 

1.08E+07 

1.07E+07 

1.08E+07 

1.07E+07 

1.07E+07 

1 .06E+07 

FCSW 

63500 

56000 

56000 

54000 

56000 

56000 

56000 

67000 

DSW 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

KDEW 

1 

1 

1 

1 

I 

1 

1 

1 

KDFW 

1 

1 

1 

1 

1 

1 

1 

1 
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Table 8. Continued 


Variable 

720 

727 

737 

747 

DC-8 

MD-11 MD-83 

L-1011 


Fuselage 


Geometric parameters 

CLBRl 

i.i 

1.1 

u 

1.1 

1.1 

1.1 

1.1 

1.1 

ICYL 

1 

1 

i 

1 

1 

1 

1 

i 

Structural concept 

CKF 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

5.24 

EC 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

2.36 

KGC 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

0.368 

KGW 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

0.505 

Material properties 

FTS(T/B) 

58500 

58500 

58500 

58500 

64000 

58500 

58500 

58500 

FCS(T/B) 

54000 

54000 

54000 

54000 

39000 

54000 

54000 

54000 

ES(T/B) 

1.07E+07 

1 .07E+07 

1 .07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1 .07E+07 

EF(T/B) 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1.07E+07 

1 .07E+07 

DS(T/B) 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

DF(T/B) 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

0.101 

TMG(T/B) 

0.04 

0.04 

0.036 

0.071 

0.05 

0.055 

0.055 

0.075 

KDE 

1 

1 

1 

1 

1 

1 

1 

1 

KDF 

1 

1 

1 

1 

1 

1 

1 

1 

Loads 

AXAC 

0 

0 

0 

0 

0 

0 

0 

0 

CMAN 

1 

1 

1 

1 

1 

1 

1 

1 

ILOAD 

3 

3 

3 

3 

3 

3 

3 

3 

PG(T/B) 

12.9 

12.9 

11.25 

13.65 

13.155 

11.5 

12.5 

12.6 

WFBUMP 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

WFLAND 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 

0.9 
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Table 8. Concluded 


Variable 

720 

727 

737 

747 

DC-8 

MD-11 MD-83 

L-1011 


Landing gear 


VSINK 

10 

10 

10 

10 

10 

10 

10 

10 

STROKE 

1.67 

1.167 

1.167 

2.21 

1.375 

1.9 

1.67 

2.17 

CLRG1 

O.I33 

0.1306 

0.145 

0.1131 

0.108 

0.141 

0.055 

0.161 

CLRG2 

0.51 

0.5896 

0.5254 

0.466 

0.499 

0.57 

0.597 

0.56 

WFGR1 

0.00389 

0.00725 

0.0052 

0.0047 

0.0311 

0.0031 

0.004 

0.005 

WFGR2 

0.036 

0.0738 

0.0382 

0.0398 

0.0742 

0.0058 

0.035 

0.044 

IGEAR 

2 

2 

2 

2 

2 

2 

2 

2 

GFRL 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

0.001 

CLRGW1 

0.1675 

0.1736 

0.1846 

0.064 

0.14 

0.2 

0.148 

0.232 

CLRGW2 

0 

0 

0 

0.1844 

0 

0 

0 

0 

Tails 

IT AIL 

1 

1 

1 

1 

1 

1 

i 

1 

Weights 

WTFF 

0.3263 

0.2625 

0.156 

0.262 

0.418 

0.336 

0.2795 

0.246 

CBUM 

1 

1 

1 

1 

I 

1 

1 

I 

CLAN 

0.813 

0.859 

0.972 

0.791 

0.7164 

0.7137 

0.9143 

0.851 


Factors 


ECHRENK 

1 

1 

1 

I 

1 

1 

1 

1 

ICOMND 

1 

1 

1 

1 

1 

1 

1 

1 

WGNO 

1 

1 

1 

I 

1 

1 

1 

1 

SLFMB 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

1.2 

WMIS 

0 

0 

0 

0 

0 

0 

0 

0 

WSUR 

0 

0 

0 

0 

0 

0 

0 

0 

WCW 

1 

1 

1 

1 

1 

1 

1 

1 

WCA 

0 

0 

0 

0 

0 

0 

0 

0 

NWING 

40 

40 

40 

40 

40 

40 

40 

40 


38 


$PDCYLIN 


PS=1„ 

EFFC=1.03, 

KDEW=1.0, 

TMGW=.02, 

ESW=10.7E06, 

KDFW=1.0, 

EFFW=.656, 

FCSW=54000„ 

DSW=0.101, 

ISTAMA=2, 

CS 1=0.088, 

CS2=0.277, 


CLAQR=.001, 

CKF=5.24, 

IFUEL=2, 

EC=2.36, 

CWMAN=1.0, 

KGC=.368, 

CF=6.25E-05, 

KGW=.505, 

FTST = 4*58500.,8*0., 
FCST = 4*54000.,8*0„ 
EST = 4*10.70E06,8*0., 
EFT = 4*10.70E06,8*0., 
DST = 4*. 101,8*0., 

DFT = 4*. 101,8*0., 
TMGT = 4*.07 1,8*0., 
KDE = 0.9, 

FTSB = 4*58500.,8*0„ 
FCSB = 4*54000., 8*0., 
ESB = 4*10.70E06,8*0., 
EFB = 4*10.70E06,8*0., 
DSB = 4*. 101,8*0., 

DFB = 4*. 101,8*0., 
TMGB = 4*.07 1,8*0., 
KDF = 0.9, 

CLBR 1=1.1, 

ll 

_l 

u 



KCONT = 12*4, 

KCONB = 12*4, 



AXAC=0., 
CMAN=1.0, 
WFBUMP=0.00 1 , 

CBUM=1.0, 

ILOAD=3, 

WFLAND=0.9, 

CLAN=0.791, 
PGB = 12*13.65, 

PGT = 12*13.65, 

WTFF=0.262, 




VSINK=10.0, 
WFGR 1=0.0047, 
CLRGW 1=0.064, 

STROKE=2.21, 
WFGR2=0.0398, 
CLRGW2 =0.1844, 

CLRG1=. 1 131, 
IGEAR=2, 

CLRG2=0.466, 

GFRL=0.001, 

ITAIL=1, 




ISCHRENK=1, 

WMIS=0„ 

NWING=40, 

ICOMND=l, 

WSUR=0„ 

WGNO=1.00, 

WCW=1.0, 

SLFMB=1.2, 

WCA=0.0, 


SEND 


Figure 16. PDCYLIN namelist for 747-21 P. 



SPAN 

BS 

ROOTC 

TIPC 

TAPER 

TRATWR 

TRATW 

TGAML 

GAMT GAMS 

VWING 

WFUEL 

DENW 


FT 

FT 

FT 

FT 




DEG. 

DEG. DEG. 

FT3 

LBS 

LB/FT3 


195.1006 

1 14.493 

41.4650 

11.0904 

0.247 

0.179 

0.078 

42.98067 30.30835 40.200 

28673.285 

186806.00 8.226 


WING 

CHORD 

LENGTH 

LENGTH 

BEND 


WEB 

COVER 

WEB CGAGE 

WGAGE 

UNITWr 

UNITWT 

NJW 

STATION 


PRIME 


MOM 


SPACE 

THICK 

THICK THICK 

THICK 

COVERS 

WEBS 


FT 

FT 

FT 

FT 

FT-LBS 


IN 

IN 

IN IN 

IN 

LB/FT2 

LB/FT2 


114.493 

11.0904 

5.4566 

5.4566 

4288. 


0.4491 

0.0543 

0.03960 0.0200 

0.02000 

0.7904 

0.5760 

3 

1 11.631 

11.8497 

5.8733 

5.8733 

27083. 


0.7541 

0.0543 

0.03960 0.0200 

0.02000 

0.7904 

0.5760 

3 

108.769 

12.6091 

6.2901 

6.2901 

77924. 


1.0254 

0.0543 

0.03960 0.0200 

0.02000 

0.7904 

0.5760 

3 

105.906 

13.3684 

6.7069 

6.7069 

162376. 


1.2780 

0.0543 

0.03960 0.0200 

0.02000 

0.7904 

0.5760 

3 

103.044 

14.1278 

7.1236 

7.1236 

284756. 


1.5198 

0.0543 

0.03960 0.0200 

0.02000 

0.7904 

0.5760 

3 

100.182 

14.8872 

7.5404 

7.5404 

448647. 


1.7549 

0.0680 

0.03960 0.0250 

0.02000 

0.9893 

0.5760 

5 

97.319 

15.6465 

7.9572 

7.9572 

657153. 


1.9857 

0.0874 

0.03960 0.0322 

0.02000 

1.2710 

0.5760 

5 

94.457 

16.4059 

8.3739 

8.3739 

912984. 


2.2139 

0.1071 

0.03960 0.0394 

0.02000 

1.5579 

0.5760 

5 

91.595 

17.1653 

8.7907 

8.7907 

1218559. 


2.4403 

0.1268 

0.03960 0.0467 

0.02000 

1.8442 

0.5760 

5 

88.732 

17.9246 

9.2075 

9.2075 

1576052. 


2.6658 

0.1462 

0.03960 0.0538 

0.02000 

2.1260 

0.5760 

5 

85.870 

18.6840 

9.6242 

9.6242 

1987387. 


2.8907 

0.1650 

0.03960 0.0607 

0.02000 

2.4002 

0.5760 

5 

83.008 

19.4434 

10.0410 

10.0410 

2454326. 


3.1155 

0.1832 

0.03960 0.0674 

. 0.02000 

2.6648 

0.5760 

5 

80.145 

20.2027 

10.4578 

10.4578 

2978392. 


3.3404 

0.2007 

0.03960 0.0738 

0.02000 

2.9183 

0.5760 

5 

77.283 

20.9621 

10.8745 

10.8745 

3561009. 


3.5655 

0.2173 

0.03960 0.08Q0 

0.02000 

3.1600 

0.5760 

5 

74.421 

21.7215 

11.2913 

11.2913 

4203376. 


3.791 1 

0.2330 

0.03960 0.0858 

0.02000 

3.3893 

0.5760 

5 

71.558 

22.4808 

11.7081 

11.7081 

4906548. 


4.0172 

0.2479 

0.03960 0.0912 

0.02000 

3.6061 

0.5760 

5 

68.696 

23.2402 

12.1248 
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Figure 17. PDCYL wing weight output for 747-21 P 
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Figure 18. PDCYL fuselage weight output for 747-21 P. 
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Figure 18. Concluded. 
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Figure 19(b). Wing weight estimation comparison for 747-21 P. 
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Appendix B - High-Altitude Study 

Description 

A study was made to estimate the wing weight of a scaled 
version of an existing propeller-driven high-altitude drone 
aircraft. This aircraft, termed the Strato7, is modeled as an 
enlarged version of the existing Perseus-a3. PDCYL was 
used to validate the wing weight estimation returned by . 
ACSYNT. 

The wing of the Strato7 incorporates a single hollow, 
cylindrical carbon-fiber/epoxy spar placed at the leading 
edge. The strength of the cover is assumed negligible. No 
fuel is carried in the wing, while propulsion and landing 
gear are mounted on the fuselage. The layout of the 
Strato7 is shown in figure 20. 


Input 

Fuselage weight estimation is not considered for the 
Strato7. An example of the ACSYNT input for the 
Strato7 wing weight estimation is shown in figure 21 . 

The corresponding PDCYLIN namelist for the case where 
the ratio of structural chord to total chord is 0.2 is shown 
in figure 22. 

Output 

Wing weight as a function of the ratio of structural chord 
to total chord is shown in figure 23. The wing weight 
estimated by ACSYNT is 789 pounds. PDCYL matches 
this wing weight when the ratio of structural chord to total 
chord is approximately 0.25. Nonoptimum weight was 
not considered in this analysis. In order to estimate 
nonoptimum weight, nonoptimum factors would need to 
be recomputed for this type of aircraft. 



44 



TRANSPORT 

4 2 2 570 570 0 0 0 2 1 7 0 
0.00010 0.6 10000.0 
12 3 6 
1 6 
1 6 

*** GEOMETRY *** 

$FUS BDMAX = 3.00, BODL =24.358, FRAB =2.01 
FRN =2.15, SFFACT = 1.082664, ITAIL = 1, 

OUTCOD = 3, SEND 

SWING AR = 23.328, AREA = 500.00, DIHED = 5.0, 

FDENWG = 0.0, LFLAPC = 0.00, SWEEP = 0.00, 

SWFACT = 1 .0, TAPER = 0.695, TCROOT = 0. 14, 

TCTIP =0.14, TFLAPC = 0.0, WFFRAC = 0.0, 

XWING =0.5664, ZROOT = 1.00, KSWEEP = 2, 

SEND 

SHTAIL AR =5.96, AREA =23.09, SWEEP =5.00, 

SWFACT = 1.0, TAPER =0.682, TCROOT = 0.08, 

TCTIP = 0.08, XHTAIL = 1 .25, ZROOT = 1 .25, 

KSWEEP = 0, SIZIT = T, HTFRAC = -0.20, 

CVHT = -2.70560, SEND 

SVTAIL AR =3.08, AREA =17.69, SWEEP =5.00, 

SWFACT = 1 .00, TAPER = 0.554, TCROOT = 0.08, 

TCTIP =0.08, VTNO = 1 .0, XVTAIL= 1.39, 

YROOT =0.00, ZROOT = 1.0, KSWEEP = 0, 

SIZIT =T, VTFRAC = -0.20, CVVT =-0.59909, 

CGM = 0.40, SEND 

SCREW NCREW =0, SEND 

SFUEL DEN =63.78, FRAC =1.00, SEND 

SFPOD 

DIAM = 2., LENGTH = 2., X = 0.592 

THETA = 90.0, SYMCOD = 1, SOD = -2, 

SEND 

SENGINEN =1, SEND 

J J 5 J J J ^ J J J J J ^ J ^ J ^ J J J ^ J ^ ^ ^ J ^ J J $$$$$$$$$ $$$$ ^ ^ J J ^ ^ $$$ $ J J J J J J J ^ J J J J 

STRDATA CRMACH = .40, QMAX = 70.45, DESLF = 2.5, ULTLF = 3.75, 

WFUEL = 392.0, WFEXT = 0.0, WFTRAP = 0.1, FRFURE = 0.0, 

IPSTO 1 = 5, TIMTOl = 0.0, IPST02 = 2, TIMT02 = 1 .0, 

IPSLND = 5, MODLND = 7, VMRGLD = 1 .2, WKLAND = 0.1, 

IBREG =0, IENDUR = 0, WCOMBP = 0.6, MMPROP = 7, 

NCODE =0, NCRUSE = 1, RANGE = 100.0, LENVEL =.FALSE., 

NLEGCL = 30, NLEGLO = 4, SEND 

2 

MACH NO. ALTITUDE HORIZONTAL NO. VIND 

PHASE START END START END DIST TIME TURN "G"S WKFUEL M IP IX W B A P 


CLIMB .414 -1 100 -1 0.0 0.0 0.0 -1 1.0000 7 41 0 0 0 0 0 

LOITER .400 -1 90000 -1 0.0 10.0 0.0 0.0 1.0000 7 21 0 0 0 0 0 


***** AERODYNAMICS ***** 

$ ACHAR ABOSB=0.074, ALMAX=20.0, AMC=12.0, ALELJ=3, ISMNDR=0, SFWF=0.99, 
SMNSWP = 0.01, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.51, 

CLOW =0.3821,0.3828,0.4266,0.4809,0.4849,0.4888,0.4946,0.5147,0.5502,0.5692, 

CMO = -. 1 59 1 1 596,-. 153 1 1 466,-. 1 502,-. 1 538,-. 158 1 ,-.1653,-. 1 749,-. 1 823, 

SEND 

SAMULT FCDW=1.1,$END 

SATRIM FVCAM =0.9183,0.9244,0.9538,0.9196,0.9230,0.9276,0.9349,09345,0.9264,0.9247, 
FLDM =1.021 1,1.0254,1.0200,1.0139,1.0200,1.0232,1.0234,1.0205,1.0226,0.8790, 

FLAP! = 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 

ITRIM = 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 

CGM=0.40, CFLAP=0.0, SPANF=0.0, IVCAM=1, ALFVC=5,0, SEND 


Figure 21. ACSYNT input for Strato/. 



$ ADET ICOD= 1 , IPLOT= 1 , N ALF= 1 0, NMDTL= 1 0, 

ALIN= -6.8, 0.0, 1.0, 2.0, 4.0, 6.0, 8.0, 10.0, 12.0, 14.0, 

ALTV = 22740.,37475.,50 1 3 1 .,6 1 224. ,7 1 097.,79992.,86 1 29. ,90000., 

SMN = 0.085, 0.1 19, 0.161, 0.210, 0.266, 0.328, 0.379, 0.400, 

ISTRS= 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

ITB= 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

ITS= 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 

SEND 

SADR AG CDBMB= 1 0*0.0, 

CDEXTR= 10*0.0, 

CDTNK= 10*0.00, 

SEND 

SATAKE DELFLD=0.0, DELFTO=0.0, DELLED=0.0, DELLTO=0.0, ALFROT=8.0, SEND 

SAPRINT KERROR=2, SEND 


Spark Ignition Internal Combustion Engine with Triple Turbocharging 
SPCONTR HNOUT = 0.,3 1001. ,50 131. ,79992., 90000., 

SMNOUT = 0.0, 0.085, 0.161, 0.328, 0.400, 

NOUTPT = 5, SEND 

SPENGIN ENGNUM = 1 , NTPENG = 4, ESZMCH = 0.00, 
ESZALT = 0., XNMAX = 7200.0, HPENG = 1 15.0, 
SWTENG = 6.0, HCRIT = 90000., FSFC = 1 .0, 

SEND 

SPROP AF =125.0, BL =2, CLI =0.5, 


DPROP = 17.88, FPRW = 0.087437, FTHR =1.0, 
NTPPRP = 12, PSZMCH = 0.00, PSZALT = 0., 

SEND 

SPGEAR GR =7.43, ETR =.95, FGRW =0.2476234, 
GRSND = 14.86, SEND 

SPENGNC XLENG = 1 .5, RLENG = 1.0, DIA1 =1.0, 

FT =0.0, FRPN =1.0, FRBT =2.0, 

NBDFT = 0.3, ANACHP = 0., DQ = 0.024, 

SEND 

TRANSPORT 

*** WEIGHTS *** 


SOPTS WGTO = 3000.0, KERROR = 2, 

SLOPE(l) =0.47970, TECHI(l) = 0.85, 
SLOPE(2) =0.97945, TECHI(2) = 0.85, 
SLOPE(3) =0.64225, TECHI(3) = 0.85, 
SLOPE(4) =0.85841, TECHI(4) = 0.85, 
SLOPE(6) =0.70145, TECHI(6) = 0.85, 
SLOPE(7) =0.85396, 

SLOPE(8) =0.55290, TECHI(8) = 0.85, 
SLOPE(9) = 1.89582, TECHI(9) = 0.85, 
SLOPE( 1 0) = 1.49618, 

SLOPE(l 1) = 0.19543, 


SLOPE( 12) = 0.48091, 
SLOPE(13) = 3.68569, 
SLOPE(l6) = 0.02254, 
SLOPE(17) = 1.0, 


KWING = 6, 
KBODY = 3, 
SEND 

SFIXW WE =757.5, 


WFEQ = 0., 
WFS = 0., 

WPL = 0., 

SEND 


Figure 21. Concluded. 



SPDCYLIN 


PS=1., TMGW=.05, EFFW=.605, 

EFFC=1.108, ESW=12.9E06, FCSW=75000., DSW=0.058, 

KDEW=1.0, KDFW=1.0, 

ISTAMA=2, CS 1=0.01, CS2=0.75, 

CLAQR=.001, IFUEL=I, CWMAN=1.0, CF=6.25E-05, 

CKF=5.24, EC=2.00, KGC=.368, KGW=.505, 

FTST = 4*58500.,8*0„ 

FCST = 4*54000.,8*0„ 

EST = 4*10.70E06,8*0., 

EFT = 4*30.0E06,8*0., 

DST = 4*. 101,8*0., 

DFT = 4*.292,8*0„ 

TMGT = 4*.03,8*0„ 

KDE=0.9, 

CLBRU1.I, 

KCONT = 12*4, KCONB = 12*4, 

AXAC=0„ CBUM=1.0, CLAN=0.93, 

CMAN=1.0, lLOAD=3, PGB = 12*1 1.5, PGT = 12*1 1.5, 

WFBUMP=0.00 1 , WFLAND=0.9, 

WTFF=0.07, 

VSINK=10.0, STROKE=1.0, CLRG1=.395, CLRG2=0.5, 

WFGR 1 =0.003 1 , WFGR2=0.0058, IGEAR=1, GFRL=0.00l, 

CLRGW1 =0.20, CLRGW2 = 0.0, 

ITAIL=1, 

ISCHRENK=1, lCOMND=l, WGNO=1.00, SLFMB=1.2, 

WMIS=0., WSUR=0„ WCW=1.0, WCA=0.0, 

NWING=40, 

SEND 

Figure 22. PDCYLIN namelist input for Strato7. 


FTSB = 4*58500.,8*0., 
FCSB = 4*54000.,8*0„ 
ESB = 4* 10.70E06,8*0„ 
EFB = 4*30.0E06,8*0., 
DSB = 4*.101,8*0„ 

DFB = 4*.292,8*0„ 
TMGB = 4* .03, 8*0., 
KDF=0.8, 

ICYL = l. 
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